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A  simplified  model  is  developed  to  analyze  the  high  frequenev 
electromagnetic  scattering  from  a  ship  at  sea.  The  ship  is  modelled 
by  an  elliptic  cylinder  to  represent  the  superstructure,  superim¬ 
posed  on  a  portion  of  an  elliptic  cone  to  represent  the  hull  struc¬ 
ture.  The  Geometrical  Theory  of  Diffraction  in  conjunction  with 
an  efficient  specular  point  tracking  algorithm  is  employed  to  cal¬ 
culate  the  scattering  from  this  ship  model.  The  sea  surface  is  — *■> 
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replaced  by  a  ""rectangular  patcir  model  which  divides  the  antenna 
illuminated  sea  surface  into  square  sub-areas.  The  scattering  from 
each  sub-area  is  represented  by  a  polarization  transforming  reflec¬ 
tion  coefficient  derived  from  the  bistatic  radar  cross-section  of 
the  patch.  .  The  cross-section  used  here  is  that  of  the  composite 
model  in  wmch  the  scattering  cross-section  of  the  ocean  is  derived 
from  the  sum  of  the  geometrical  optics  cross-section  (which  pre¬ 
dominate;/  in  the  forward  scatter  direction  and  is  controlled  by 
the  laraje  scale  slope  distribution  of  the  ocean)  and  the  slope  av¬ 
eraged  Bragg  scatter  (which  predominates  in  the  back-scatter  di¬ 
rection  and  is  controlled  by  the  height  spectrum  of  the  ocean). 

Mn  general,  there  are  three  basic  scattering  mechanisms  that 
contribute  to  the  total  signature  as  the  radar  system  traverses 
the  ship  target.  These  are  the  direct  back-scattered  signal  from 
the  sea  surface  (the  so-called  clutter  component),  the  ship-sea 
double-bounce  interaction  term,  and  the  signature  scattered  by  the 
ship  alone.  A  number  of  signatures,  for  various  trajectories,  an¬ 
tenna  patterns,  and  polarizations  have  been  calculated,  and  are 
used  to  illustrate  the  effect  of  the  aspect  angle  of  the  trajectory, 
the  curvature  of  the  hull,  etc.  on  the  detectability  of  the  ship 
in  clutter,  and  the  relations  between  the  position  of  the  ship  and 
the  position  of  the  peak  signature  voltage.  __ 
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CHAPTER  I 
INTRODUCTION 


When  an  antenna  is  transmitting  or  receiving  in  the  vicinity 
of  the  sea  surface,  the  desired  signal  is  usually  contaminated  by 
the  presence  of  signal  components  that  have  been  scattered  from 
the  sea.  This  form  of  Interference  generally  degrades  system  per¬ 
formance  unless  the  system  can  be  designed  to  take  advantage  of 
any  differences  between  the  properties  of  the  sea-scattered  and 
the  desired  signals.  The  desired  signal,  in  this  study.  Is  the 
radar  return  due  to  the  presence  of  a  ship.  In  order  to  determine 
the  characteristics  of  this  signal  and  to  obtain  a  better  under¬ 
standing  of  the  sea  and  ship  scattering  mechanisms,  we  have  chosen 
to  study  the  scattered  signals  through  the  use  of  a  mathematical 
model.  Such  a  model  should  be  complex  enouqh  to  realistically  re¬ 
produce  all  the  significant  features  of  the  real  situation,  while 
remaining  computationally  tractahle.  The  basic  form  of  the  model 
developed  here  is  illustrated  in  Fiqure  1  where  the  ship  Is  modelled 
by  a  pair  of  elliptic  cylinders  which  are  taken  to  represent  the 
ship's  hull  and  superstructure.  Although  this  structure  is  rela¬ 
tively  simple,  the  model  has  the  potential  of  handlinq  considerably 
more  complex  structural  features  without  any  fundamental  change 
in  the  computational  process.  Also  in  Figure  1  the  actual  sea  surface 
is  replaced,  in  +he  model,  by  a  set  of  square  patches  which  divides 
area  illuminated  by  the  antenna  into  sub-areas. 


Figure  1.  Basic  ship-sea  model. 


As  the  radar  system  traverses  the  ship  target,  a  signature 
(typically  the  detector  voltage  versus  time)  Is  generated.  There 
are  generally  three  principal  contributions  to  the  total  signature. 
The  first  of  these  Is  the  direct  backseat tered  signal  from  the  sea 
surface,  usually  referred  to  as  the  sea  clutter  signal.  The  second 
is  the  "double-bounce"  corner  term  formed  by  any  scattering  path 
from  transmitter  to  sea  surface  to  ship  and  back  to  receiver  (and 
vice  versa).  The  third  type  of  contribution  Is  the  scattered  signal 
from  the  ship  alone,  predominantly  due  to  reflection  and  diffraction 
from  the  ship's  surfaces  and  edges.  The  object  of  the  work  reported 
here  Is  the  prediction  of  these  three  principal  signatures.  They 
will  be  analyzed  by  two  separate  techniques,  namely  the  statisti¬ 
cally  rough  surface  theories  for  the  sea  scattering,  and  the  Geo¬ 
metrical  Theory  of  Diffraction  (GTD)  for  the  ship  scattering. 

The  calculation  of  scattering  of  the  electromagnetic  field 
hv  the  sea  surface,  which  Includes  hoth  forward  and  back  scattering, 
is  itself  based  on  a  model.  In  this  case  the  "composite"  model  , 
in  which  a  slightly  rough  surface  Is  superimposed  on  a  larger  scale 
structure.  Perturbation  theory  Is  used  to  find  the  Bragg  scatter 
from  the  small  scale  structure,  and  the  methods  of  Geometrical  Optics 
are  employed  for  scattering  from  the  large  scale  ( gravity  wave) 
structure.  The  scattered  fields  are  then  formulated  In  terms  of 
the  bistatic  radar  cross  section  of  each  patch,  as  computed  by  the 
above  two  statistical  rough  surface  theories.  Chapter  IV  discusses 
the  characteristics  of  the  sea  surface  scattering  In  more  detail. 

Both  the  Physical  Optics  method  and  the  Geometrical  Theory 
of  Diffraction  (GTD)  have  been  Investigated  for  the  calculation 
of  the  electromagnetic  scattering  from  the  ship.  It  was  found  that 
the  GTD  is  not  only  more  effective  in  computing  the  scattering  from 
such  objects  as  elliptic  cylinders,  plates,  rectangular  blocks, 
etc.,  but  also  allows  one  to  gain  more  physical  Insight  Into  the 
various  scattering  and  diffraction  mechanisms  Involved.  Therefore 
the  method  of  GTD,  which  Is  described  In  Chapter  III,  was  chosen 
to  formulate  the  problem  of  scattering  from  the  ship. 

Finally,  in  order  to  determine  the  double-bounce  ship-sea  Inter¬ 
action  term,  we  have  again  used  the  GTD,  In  which  the  sea-scattered 
fields  from  each  patch  are  represented  as  a  new  source.  To  carry 
out  this  procedure.  It  is  essential  to  locate  accurately  the  re¬ 
flection  point  on  the  surface  of  the  ship's  hull.  Since  huge  num¬ 
bers  of  patches  might  be  Involved  If  the  area  Illuminated  bv  the  an¬ 
tenna  Is  large,  a  great  deal  of  time  could  be  spent  in  just  com¬ 
puting  the  reflection  points.  For  this  reason,  a  verv  efficient 
numerical  search  technlgue  is  developed  to  track  the  reflection 
point  as  the  source  moves  from  one  sub-area  to  another.  Chapter 
V  has  a  detailed  discussion  of  this  technlgue. 


The  combined  signatures  for  the  complete  model  are  considered 
in  Chapter  VII.  In  this  chapter  the  hull  form  is  changed  from  an 
elliptic  cylinder  to  an  elliptic  cone  so  that  the  geometrical  shape 
of  the  ship  is  physically  more  realistic.  Also  in  the  same  chapter, 
the  antenna  pattern  is  modified  from  a  pencil  beam  to  a  conical 
beam,  so  the  effects  of  different  antenna  patterns  on  the  signatures 
can  be  analyzed.  In  the  concluding  chapter,  the  general  problem 
of  detectability  of  a  ship  target  in  clutter,  and  some  directions 
for  further  study  are  considered. 
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CHAPTER  II 

THE  THREE  DIMENSIONAL  SIMPLIFIED  SHIP-SEA  MODEL 


It  is  the  purpose  of  this  chapter  to  give  a  basic  geometrical 
description  of  the  simplified  ship-sea  model.  In  this  model,  the 
radar  system  is  represented  by  a  transmit/receive  antenna  of  speci¬ 
fied  pattern,  moving  along  a  specified  straight-line  trajectory. 
Because  the  interest  is  in  radar  wavelengths  much  smaller  than  the 
characteristic  dimensions  of  the  ship,  and  in  trajectories  which 
transmit  at  relatively  short  ranqe,  the  problem  is  formulated  here 
in  terms  of  the  Geometrical  Theory  of  Diffraction.  The  ship,  as 
shown  in  Figure  ?,  is  first  modelled  by  a  single  elliptic  cylinder 


Figure  2.  Simplified  models  of  the  ship's  hull 
and  superstructure. 
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representing  the  hull.  Later  on  a  smaller  elliptic  cylinder  is 
added  to  the  top  to  simulate  the  superstructure;  and,  finally,  the 
ship's  hull  is  modified  to  have  the  shape  of  an  elliptic  cone.  Since 
the  elliptic  cylinders  are  three  dimensional  structures,  a  3-D  anal¬ 
ysis  is  necessary. 

The  basic  geometry  for  the  ship-sea  interaction  calculation 
is  shown  in  Figure  3,  in  which  the  antenna  is  traveling  in  an  ar¬ 
bitrary  direction/!  with  its  main  beam  axis  (bore- sight)  pointed 
in  the  direction’s  and  is  illuminating  the  sea  surface  over  an  ef¬ 
fective  area  of  WxL.  This  effective  area,  which  is  dependent  on 
the  transmitting  antenna  pattern,  is  broken  up  into  smaller  sub- 
areas.  The  scattering  from  each  sub-area  is  entered  into  the  GTD 
program  as  a  polarization  dependent  scattering  matrix  element  with 
scattering  coefficients  referred  to  the  center  of  the  sub-area. 

Each  sub-area  generates  a  direct  backscattered  field  to  form  the 
sea  clutter  signal.  The  forward  scattered  field  which  is  also  gen¬ 
erated  from  each  sub-area  then  appears  to  come  from  a  new  source 
which  in  turn  illuminates  the  metallic  cylinder,  and  from  which 
it  is  reflected  to  the  receiver.  The  form  of  this  doubly  reflected 
field  (the  ship-sea  interaction  term)  constitutes  one  of  the  major 
objectives  of  this  research,  because  it  represents  one  of  the  im¬ 
portant  contributions  to  the  overall  target  signature.  Even  though 
this  component  of  the  scattered  field  is  somewhat  smaller  than  the 
signal  scattered  by  the  ship  alone,  it  is,  in  many  instances,  a 
stronger  signal  than  that  due  to  the  sea  clutter.  Since  this  doubly 
reflected  field,  whose  existence  is  sometimes  overlooked,  occurs 
earlier  than  the  ship  scattered  signal,  it  can  play  an  important 
role  in  the  overall  system  design,  particularly  when  accurate  range 
estimates  are  required. 

The  antenna  radiation  pattern,  as  will  become  apparent  later, 

can  also  play  an  imoortant  role  in  the  overall  signature  of  the 

ship-sea  model.  However,  for  convenience,  the  antenna  pattern  is 
first  chosen  to  be  a  pencil  beam  radiated  from  a  rectangular  aper¬ 
ture  with  Gaussian  amplitude  distribution  in  both  planes.  The  rec¬ 
tangular  aperture  with  its  coordinate  system  and  the  radiated  fields 
are  given  in  Appendix  A.  Although  a  near  field  analysis  is  required 
to  specify  the  range  dependence  of  the  radiated  field  (because  the 

antenna  may  be  in  the  near  zone  of  the  ship  at  the  range  of  inter¬ 

est),  it  can  be  shown  that  the  angular  dependence  of  the  antenna 
pattern  is  equivalent  to  a  far  field  pattern. 

Similarly,  the  details  of  the  system  electronics  (detector 
type,  range  gate,  doppler  filters,  etc.)  can  play  an  important  role 
in  the  detectability  of  targets  in  clutter.  For  the  purposes  of 
this  study,  however,  we  have  assumed  that  the  emitted  signal  Is 
a  c.w.  waveform  proportional  to  [Re  V  exp(2njf  t)]  where  V  Is 
a  transmitter  reference  voltage,  f  is  the  carrier  frequency,  and 
Re  designates  the  real  part  of  theccomplex  expression.  The  model 
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(b)  PERSPECTIVE  SIOE  VIEW 


Figure  3 


Geometry  of  the  cylinder  and  rectangular 
patch  model. 
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represents  the  received  signal  as  a  sum  of  contributions  with  well 
defined  total  path  length  r  (e.g.,  for  double  bounce  path,  r  Is 
the  round  trip  path  length  from  transmitter,  to  patch  center, nto 
reflection  point  back  to  the  recelverl.  Thus  the  received  signal 
has  the  form 

Re  l  Vn  expf?irffct  -  jkrp  +  6nl 
n 

where  k  =  ?n/X,  V  Is  an  amplitude  factor,  and  is  a  phase  (pos¬ 
sibly  random!  which  depend  on  the  scattering  process  for  the  par¬ 
ticular  contribution.  The  detector  is  assumed  to  be  of  the  homodyne 
type  with  low  pass  filter,  so  that  the  detected  signal  V.(t)  has 
the  form  0 

Vd(t)  =  Re  l  Vn  exp(-jkrn+6n). 

This  is  the  signal  which  is  computed  in  the  subsequent  sections. 

Note  that  the  doppler  shift  and  the  general  time  dependence  of  the 
waveform  are  contained  implicitly  in  the  constantly  changing  path 
lengths  r  . 
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CHAPTER  III 

INTRODUCTION  TO  THE  GEOMETRICAL  THEORY  OF  DIFFRACTION  (GTD) 


In  this  chapter  the  basic  GTD  solutions  used  to  analyze  the 
simplified  ship  model  are  reviewed.  The  GTD  is  a  high  frequency 
technique  that  allows  a  complicated  structure  to  be  approximated 
by  basic  shapes  such  as  flat  plates,  curved  wedges  and  convex  curved 
surfaces.  It  is  a  ray  optical  technique  and  therefore  allows  one 
to  qain  some  physical  insight  Into  the  various  scattering  and  dif¬ 
fraction  mechanisms  Involved.  Consequently,  one  is  able  to  quickly 
seek  out  the  dominant  or  significant  scattering  mechanism  for  a  given 
qeometrical  configuration.  Several  fundamental  scattering  mech¬ 
anisms  are  used  in  the  GTD  calculation  for  the  simplified  ship  model. 
They  are:  1.  Reflected  field  from  a  curved  surface,  ?.  Diffrac¬ 
tion  hy  a  straight  wedge,  and  3.  Diffraction  by  a  curved  wedge. 

These  scattering  mechanisms  are  individually  discussed  below. 

A.  Reflected  Field  from  a  Curved  Surface 


The  scattering  problem  studied  here  is  that  of  the  reflection  of 
a  given  electromagnetic  field  E ' ( QR )  from  a  curved  surface,  as  shown 
in  Figure  A.  The  reflected  electric  field  Er(s)  at  the  observation 
point  is  given,  in  the  geometrical  optics  approximation,  by? 


Figure  4.  Reflection  from  a  curved  surface. 


where  F  fQ.,1  is  the  incident  electric  field  at  the  point  of  reflec¬ 
tion  0B  due  to  a  point  source  at  P.  the  reflection  point  QR  Is 
found  from  the  law  of  reflection,  which  requires  that  the  angle  of 
incidence  he  equal  to  the  anqle  of  reflection.  It  is  determined 
hy  the  relation 

a  \  >  \ 

-1  •  n  *  s  •  n  (?al 

N  \ 

where  I  i^  the  incident  unit  vector,  s  is  the  reflection  unit  vec¬ 
tor,  and  n  is  the  surface  normal  unit  vector.  It  is  also  required 
that  all  those  unit  vectors  must  lie  in  the  same  plane  at  the  re¬ 
flection  point,  so  that 


I  x  n  »  s  x  n.  (?bl 


The  dyadic  reflection  coefficient  (T  for  a  perfect  conducting  surface 
can  he  expressed  as 


'1  'r 
e  e 


e  e 

1  A. 


(3) 


'  i 

where  e„  and  e„  are  the  incident  and  reflected  unit  vector  parallel 
to  the  plane  of  incidence  M.e.,  the  plane  containing  the  incident 
vector  and  the  surface  normal  1,  and  e*  is  the  unit  vector  per¬ 
pendicular  to  the  plane  of  incidence.  The  quantities  PT  and 
are  the  principal  radii  of  curvature  of  $he  reflected  wavefront 
at  the  reflection  point  Qp.  Kouvoum.lian  discusses  in  detail  how 
to  find  these  values  for  an  arbitrary  wavefront.  For  a  spherical 
incident  wavefront,  the  radi i  are  given  hv 


cos  p 


sin  0., 


1 

'y  o 

sin'  O,  sin'  O. 

1  * 

v 

J  COS  0 

-T,  fST 

M) 


where 


s’  is  the  radius  of  curvature  of  the  incident  wavefront  at 

Qp. 

0*  is  the  anqle  of  incidence  with  cosO1 *n*s*-n*f, 

0.  is  the  angle  between  the  incident  ray  and  the  principal 
surface  tangent  vector  e^, 

0«  is  the  angle  between  the  incident  rav  and  the  principal 
surface  tanqent  vector  e0, 

i 


o 


R,  and  R?  are  the  principal  radii  of  curvature  of  the  surface 

it  QR , 

s  is  the  distance  from  QR  to  the  observation  point  0. 
li.  Diffraction  by  a  Straight  Wedge 

Even  thouqh  the  diffraction  by  a  straight  wedge  has  no  direct 
application  to  the  ship  model  we  have  chosen,  it  is  discussed  here 
because  it  is  one  of  the  most  fundamental  diffraction  mechanisms 
in  the  GTD  analysis  and  it  is  essential  to  understand  it  before 
considering  the  diffraction  by  a  curved  wedge  in  the  next  section. 

An  asymptotic  solution  for  the  diffraction^from  a  conducting 
wedqe  was  first  obtained  by  Sommerfeld  .  Pauli4  introduced  the 
VR  function  as  a  convenient  formulation  of  the  diffraction  solution 
for  a  conducting  wedqe  of  finite  anqle.,.  Later  a  plane  wave  dif¬ 
fraction  coefficient  was  used  by  Keller  in  his  original  develop¬ 
ment.  of  the  GTD.  Recently,  however,  Hutchins  and  Kouyoumjian^  have 
presented  a  formulation  for  the  diffracted  field  which  is  signifi¬ 
cantly  more  accurate  than  the  Pauli  or  Keller  forms  in  the  transi¬ 
tion  regions  f near  the  incident  and  reflection  shadow  boundaries!. 


The  three  dimensional  wedge  diffraction  problem  is  depicted 
in  Fiqure  5.  The  incident  field  F'(s)  is  radiated  from  the  source 
located  at  the  point  s 1  (p1  ,<j> 1  ,z’ }.  This  can  be  an  arbitrary  elec¬ 
tric  or  magnetic  source  producing  a  plane,  a  cylindrical  or  a  spheri 
cal  wave  incident  on  the  wedge.  Kouyoumjian  and  Pathak'  have  re¬ 
cently  shown  that  the  diffracted  fields  may  be  written  compactly 
in  terms  of  a  dyadic  diffraction  coefficient  if  the  fields  are  ex¬ 
pressed  in  a  ray-fixed  coordinate  system.  The  ray-fixed  coordi¬ 
nate  system  is  centered  at  the  point  of  diffraction  which  is 
the  so-called  stationary  phase  point  (a  unique  point  Tor  a  given 
source  and  observation  point).  The  incident  ray  diffracts  as  a 
cone  of  rays  such  that  PQ  =  (see  Figure  5)  or 

A  A  A  A 

-I  •  2  =  0  •  z  =  cos  gQ.  (5a) 

The  relationships  between  Jthe^ortfcogon^l  unit  vectors  associated 
with  these  coordinates  (I.g’.J';  6,8  ,$)  are  shown  in  Figure  5, 
and  are  given  by  c 


I  =  &'0  x 
D*=  x  i 


(5b) 


where  S'  and  $  are  the  unit  vectors  at  the  source  and  observation 
point,  respectively. 


J 
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Figure  5.  Geometry  for  three  dimensional  wedge  diffraction. 
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The  diffracted  field  is  now  given  by 

Fd(sl  r  F1 (QeT  .  6  A(s)  e"jks  (6) 

where  s  is  the  distance  between  the  observation  point  and  diffrac¬ 
tion  point,  and?)  is  the  dyadic  diffraction  coefficient,  given  by 

^  =  'BoBons  '  Dh 


wi  th 


'’nj^nk  sinp 


cot  F  [kLa  ( p  1]  +  cot  ^  J  F  [kLa  ( p  1] 

-f  S 

cot  ('Tn")  F  [kLa’(S+)] 


(7) 


Here  the  upper  sign  (-1  corresponds  to  the  soft  boundary  condition 
(si  and  the  lower  sign  (+)  corresponds  to  the  hard  boundary  con¬ 
dition  (hi.  The  angle  variables  are  given  by  B*  =  $+0' ,  where  the 
terms  are  associated  with  the  incident  field  and  the  0+0'  terms 
are  associated  with  the  reflected  field.  The  distance  parameter 
l  will  be  defined  later.  The  function  F  [x]  is  called  the  transition 
function  and  is  given  by 


F[«]  •  2 j|j7|eJX  /  e'-”  dr 

|v/*l 

The  above  diffraction  coefficient  can  be  shown  to  be  related  to 
the  Vg  function  for  a  straight  wedge  by 

+  FT 

»  PVL.fO  5  v8a,m]J^- 

where 
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(PI 


v8a,pi  =  +  i+7/l,pi 


wi  th 


-  .i  f  k  L  +  n/4 )  . 

1  =  e'Jri  jTn - >ia  C0t(7r)e 


...  <*>  .  7 

jUa  f  0-j,  dt 


/kU 

+  [higher  order  terms] 
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The  higher  order  terms  are  negligible  for  large  kL.  Here  n  is 
defined  from  the  wedge  angle  WA  =  (2-n)ir;  also 


a  =  a1  =  1  +cosf  R-?nuN ±)  and 

is  that  positive  or  negative  inteqer  for  zero),  which  most  nearly 
satisfies  the  equations 

9nnN"  -8  =  -it  for  I 

*  IT 

?nnN+  -R  =  +n  for  I.  . 

+  TT 

For  a  straigh^  wedge  the  quantity  L  is  a  distance  parameter  given 
in  qeneral  by 

s(p’+s)p1Vsin?8 

L  =  -.■e' -t  - -  (10) 

Pg(pJ+s)(p2+s) 

where  pj,  p^  are  the  principal  radii  of  curvature  of  the  incident 
wavefront  at  the  point  of  diffraction,  and  p^  is  the  radius  of  curva¬ 
ture  of  the  incident  wavefront  in  the  plane  containing  the  incident 
ray  and  the  edge.  In  matrix  notation,  the  diffraction  coefficient 
can  be  written  as 


Ej(s) 

E*(s) 


-Ds  0  1  F  E,!(QE) 

0  -DhJ  Ul(QE) 


A(s)e" 


The  D  coefficient  fminus  siqn  between  VR  terms)  applies  to  the 
F-fieTd  component  parallel  to  the  edge  and  is  derived  from  the  can¬ 
onical  problem  with  the  acoustically  soft  boundary  condition 

fl|  .  )  =  n. 

wedge 

The  0.  coefficient  (plus  siqn  between  VR  terms)  applies  to  the  T- 
f ieldncomponent  perpendicular  to  the  edge,  and  is  derived  from  the 
canonical  problem  with  the  acoustically  hard  boundary  condition 


'  wedged 

The  quantity  A(s)  is  the  ray  divergence  factor  given  in  general 
by? 


A(s)  = 
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where  for  a  stralqht  wedqe,  p=p’  is  the  radius  of  curvature  of  the 
incident  wavefront  in  the  planecontaininq  the  incident  ray  and 
the  edqe.  For  some  specific  problems  the  diverqence  factor  is  given 

plane,  cylindrical  and  conical 
wave  incident 


spherical  wave  incident 

where  s'  is  the  distance  from  the  source  to  the  diffraction  point 
and  s  is  the  distance  from  the  diffraction  point  to  the  observ¬ 
ation  point. 

C.  Diffraction  by  a  Curved  Wedqe 

In  this  report  the  hull  of  the  ship  is  to  be  modelled  first 
by  a  finite  elliptic  cylinder.  Since  the  top  and  the  sides  of  the 
elliptic  cylinder  form  a  curved  wedge  at  their  junction,  this  con¬ 
figuration  must  also  be  considered.  The  curved  wedge  problem,  which 
is  illustrated  in  Figure  6,  is  again  analyzed  using  the  GTD  tech¬ 
nique  developed  by  Kouyoumjian  and  Pathak  .  Again,  the  diffracted 
field  from  the  curved  wedge  may  be  written  in  the  form 


Ed(s)  =  ^(Qe) 


a?) 


The  parameter  p  is  the  distance  between  the  caustic*  at  the  edg^ 
and  the  second  caustic  of  the  diffracted  ray.  This  is  given  by 


P 


•  H-s) 
-P-.  - 
sin  80 


M?1 


where  3  ,  n  ,  a  are  defined  below.  The  diffraction  coefficients 
for  thecurved  wedqe  are  extended  from  those  in  Equation  (7)  to 
allow  the  diffracted  field  to  be  continuous  at  the  incident  and 
reflected  shadow  boundaries.  This  is  accomplished  by  finding  the 
appropriate  distance  parameter  L  in  each  of  the  transition  functions 
that  makes  the  total  field  continuous.  The  diffraction  coefficient 
for  the  curved  wedge  is  given  by 


*A  caustic  is  a  location  where  the  rays  converge  and  the  geometrical 
optics  field  becomes  infinite.  It  can  be  a  point,  a  line  or  a  surface. 


-jn74 


/ 

’  sin  £  FfkL’af^')] 

)  (*t*,;h  1  =  - - A - - 

^*«o  [co5(l). 


cos(f) 


cot(1^)F[kLrV(Bt)]  ♦  cot(2^)F  [kLroafB  +  )]J  ■ 

in  which  a(B)  =  ?co$^(p/2)  and  a+(p)  =  2  cos^«n~^i 
parameter  associated  with  the  incident  field  is^given  by 


(14) 

The  distance 


,  i ,  >  i  i  .  2„ 
s(pe+s)pip?  sin  PQ 

Pe(p{+sHpp+s) 


(lba) 


ro 

The  distance  parameters  L  associated  with  the  reflecting  surface 
denoted  by  0,  or  Lrn  associated  with  the  reflecting  surface  denoted 
by  n,  are  given  hy  (see  Figure  7) 


sfPp+slpljppSin^o 
Pp(p  [+sl  (pp+s) 


(15b) 


The  parameters  p,1  andp^  are  the  principal  radii  of  curvature  of 
the  incident  wavefront  St.  the  diffraction  point  QF.  The  parameter 
p’  is  the  radius  of  curvature  of  the  incident  wavefront  at  CL  taken 
Jn  the  plane  containing  the  incident  ray  and  the  tangent  unit  vector 
e  at  the  edge  (see  Figure  6).  In  the  case  of  plane,  cylindrical 
and  conical  waves,  p^  is  infinite  and  for  the  case  of  spherical 
waves  p  =s ' .  The  parameters  pT  and  pp  are  the  principal  radii  of 
curvature  of  the  reflected  wavefront  at  (L  and  may  be  found  from 
Filiation  (4)  for  spherical  wave  incidence.  Note  that  plj,  p£,  and 
P  ^  must  be  determined  separately  for  the  o  and  n  surfaces,  in  that 
eSch  surface  has  its  own  reflection  properties.  The  parameter  p 
is  the  radius  of  curvature  of  the  reflected  wavefront  at  (L  takefi 
in  the  plane  containing  the  reflected  ray  and  the  tangent  unit  vec¬ 
tor  e  and  is  given  by 
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>(n 


ne)0 


n) 


ae  sin?Bo 


(16) 


whore  n  is  the  unit  normal  to  the  surface  either  o  or  n.  The  unit 
normal  to  the  edge  (n  )  is  directed  away  from  the  center  of  edge 
curvature,  and  a  >  Oeis  the  radius  of  curvature  of  the  edge  at 

n  ^ 
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OBSERVATION 

POINT 


(Q) 


(b)  PLANE  PERPENDICULAR  TO  •  AT  Qc 


Figure  6.  Geometry  for  three-dimensional  curved 
wedge  diffraction  problem. 


These  three  standard  GTD  solutions  will  form  the  basis  of  our 
calculations  for  the  ship  and  ship-sea  scattering  contributions 
to  the  signatures  detected  by  the  radar. 


SOURCE 


Figure  7.  Definition  of  surfaces  for  wedge 
diffraction  problem. 


CHAPTER  IV 

SEA  SCATTERING  CHARACTER ISTICS 


The  scatterinq  caused  hv  the  presence  of  the  ocean  surface 
should  he  calculated,  in  principle,  from  the  actual  state  of  the 
surface,  i.e.,  from  a  specification  of  the  heiqht  of  the  surface 
as  a  function  of  position  and  time.  This  would  he  both  difficult 
and  expensive.  At  the  present  time  there  is  not  available  a  suf¬ 
ficiently  detailed  description  of  the  ocean  surface  to  correctly 
formulate  the  microwave  scatterinq  problem.  If  there  were  such 
a  description,  the  scattering  problem  still  could  not  he  solved 
exactly  with  present  techniques.  Finally,  even  such  approximate 
techniques  as  physical  optics  would  require  a  double  inteqration 
over  the  extremely  large  surface  area  which  can  interact  with  the 
ship  and  the  radar  system.  Thus  in  the  ship- sea  model  developed 
here  the  actual  scatterinq  properties  of  the  surface  have  been  re¬ 
placed  by  a  statistically  averaged  bistatic  radar  cross-section  of 
the  ocean.  In  this  chapter  we  discuss  how  the  radar  cross-section 
approach  is  implemented  in  the  model,  and  what  particular  form  of 
radar  cross-sect  ion  is  used. 

When  usinq  the  radar  cross-section  to  characterize  the  ocean 
scatter,  care  must  he  taken  to  duplicate,  as  far  as  possible,  the 
statistical  properties  of  the  actual  received  signal.  We  have  ap¬ 
proached  this  problem  bv  dividinq  up  the  total  surface  in  the  vicinity 
of  the  ship  model  into  rectanqular  patches  (see  Fiqure  3).  Now 
the  scatterinq  from  each  patch  can  be  found  from  the  radar  cross- 
section  per  unit  area  of  the  patch,  and  its  area.  The  size  of  the 
patches  must  be  chosen  small  enough  so  that  the  angular  dependence 
of  the  cross-sect  ion  does  not  change  great  ly  from  patch  to  patch, 
but  large  enouqh  so  that  the  scattering  from  patch  to  patch  is  sta¬ 
tistically  independent.  Because  of  the  random  character  of  ocean 
surface,  the  positions  of  the  specular  points  which  contribute  to 
the  forward  cross-section  are  themselves  random.  Thus  the  minimum 
patch  size  should  contain  enough  specular  points  (say  more  than 
10)  so  that  the  signal  from  each  patch  has  a  statistically  randan 
character.  We  have  not  addressed  the  problem  of  determining  the 
optimum  patch  size  in  this  report,  but  believe,  based  on  the  con¬ 
sideration  suqgested  above,  that  the  dimension  of  the  patch  should 
be  several  times  the  wavelength  of  the  significant  ocean  wave  size, 
and  that,  the  angle  subtended  by  the  centers  of  adjacent  patches 
should  probably  not  be  much  more  than  10  degrees. 

To  specify  more  formally  the  contribution  from  a  single  patch, 
consider  the  geometry  of  Figure  A.  If  a  plane  wave  with  electric 
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field  E  exp(-jkx  sin0.+jkz  cosO.)  is  incident  on  the  patch  LxL 
with  angle  of  incidence  0^,  then^he  electric  field  in  the  scatter¬ 
ing  direction  6S,  4>$  at  range  R$  will  be  given,  in  our  model,  by 


Es  = 


IU 


*Ves**s> 


°o  E  'JkRs+J« 
~W  e 


(17) 


where  1^(0. ,0  ,4>  )  is  a  unit  vector  which  specifies  the  polarization 
state  of  thesobierved  field,  5  is  a  random  phase  angle  with  -n  < 
fi  <  it,  and  a  (0.,0  ,<j>  )  is  the  bistatic  radar  cross-section  per 
unit  area  of°theocean:  We  consider  next  the  specification  of  the 
quantity  a  . 


Figure  8.  Surface  arrangement  and  scattering  geometry. 
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The  general  problem  of  scattering  of  electromagnetic  waves 
by  the  ocean  is  an  extremely  difficult  one.  The  forces  of  wind, 
gravity,  surface  tension,  acceleration,  etc.  produce  a  spectrum  of 
ocean  waves  of  various  heights,  lengths  and  directions  of  motion 
which  collectively  behave  as  a  random  rough  surface.  The  scattering 
from  statistically  rough  surfaces  has  been  studied  extensively  in 
recent  years  and  two  important  special  cases  of  rough  surface  scat¬ 
tering,  distinguished  by  their  relative  roughness  scales,  and  by 
the  analytical  methods  used  to  obtain  the  scattered  field,  have 
particular  relevance  to  ocean  scatter.  The  currently  accepted  view 
of  the  scattering  behavior  of  the  ocean  is  contained  in  the  so-called 
"composite"  model.  In  this  model,  the  scattering  near  the  specular 
direction  M.e.,  near  normal  incidence  for  back-scatter)  is  con¬ 
trolled  by  the  slope  distribution  of  the  large  scale  structure  of 
the  surface.  This  part  of  the  scattered  field  is  calculated  by 
physical  or  geometrical  optics,  and  provides  an  explanation  not 
onlv  for  the  near-normal  incidence  back-scattering  cross-section 
but  also  for  forward  scatter  (multipath)  effects,  and  for  the  de¬ 
pendence  of  brightness  temperature  on  wind-speed  and  polarization. 

For  the  back-scatter  cross-section  away  from  normal  Incidence,  the 
scattering  mechanism  is  the  "Bragg"  scatter  from  capillary  and  short 
wavelength  components  of  the  surface.  This  contribution  is  calcu¬ 
lated  by  perturbation  theory,  and  when  utilized  with  the  known  spec¬ 
trum  of  the  ocean  it  explains  the  angular  and  polarization  depend¬ 
ence  of  the  radar  cross-section,  and  the  weak  dependence  of  this 
cross-section  on  wind  speed  and  electromagnetic  wavelength.  In 
this  composite  model,  these  two  contributions  to  the  scattered  power 
are  considered  to  be  statistically  independent,  so  that  the  total 
cross-section  may  be  written  as  the  sum, 

ao(°i’°s^s)  =  °os(9i’es’<,s)  +  ^op^i’W*  (18) 

where  oos  is  the  specular  component,  and  <aOD>  is  the  perturbation 
component  averaged  over  the  slope  distribution  of  the  large  scale 
surface  structure.  These  two  components  are  considered  separately 
as  follows. 

A.  The  Bragg  Component.;  Slightly  Rough 

Surface  Scattering  ” 

The  slightly  rough  surface  is  one  having  a  scale  of  roughness 
whose  r.m.s.  roughness  height,  h,  is  much  less  than  the  electrical 
wavelength  and  whose  surface  slopes  are  small  compared  to  unity. 

This  slightly  rough  surface  model,  based  ^n  a  perturbation  calcu¬ 
lation,  was  originally  formulated  by  R1ce?  in  1951  and  was  adapted 
to  problems  of  radar  scattering  a  few  years  later  by  Peake.  The 
perturbation  or  smallness  parameters  are  the  surface  height  -K  (x,y) 
and  the  surface  slopes  -  F>x  and  C  By  expanding  the  surface  height 
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£(x,y)  in  a  Fourier  Series,  and  applying  the  perturbation  theory 
with  the  appropriate  boundary  conditions,  one  obtains'  the  average 
scattering  cross  section  per  unit  surface  area  for  bistatic  scat¬ 
tering: 

one(O1’0s‘^s)  =  16llk4  cos20i  cos20slctncl?  S(kx*ky)  (19) 

-  where  n  and  e  are  the  scattered  and  incident  polarization  states, 

respectively, 

-  k  is  the  free  space  wave  number  (2tt/A)  of  the  electromagnetic  wave. 

-  The  angles  are  illustrated  in  Figure  8.  The  incident 

plane  is  assQmea  to  be  aligned  with  the  x-axis.  For  backscat- 
tering  fli  =  0$,  <J>S  =  it. 

-ct  is  the  scattering  matrix  element  for  a  homogeneous  dielectric 
nc  surface  and  is  a  function  of  the  polarization  states  and  material 
properties  ( see  Appendix  Bl. 

-  The  function  Sf k  ,k  )  is  physically  the  surface  height  spectral 

density,  y 

-  k  ,  k  are  the  (mechanical)  radian  wave  numbers  in  the  x  and  y 

x  directions,  respectively. 

k  =  k(sin6  cos<|>  -Sln8,),  k  =  k  sine  sin<j> 

X  5  5  *  J  5  S 

The  density  S(k  k  )  is  defined  as 
x  y 

1  ®  -JkyTy  "Jkx/Tw 

S(k  ,k  )  = - *  J7  <c(x,y)^(x',y')>  e  e  y  y  dT  dx 

y  (2lf)^  y 


where  t  =  x-x'  and  r  =  y-y', 

*  y 


s'kx*V  -77 ^//R,vVe 


e  dTx  rfTy 


where  r^tx,tv1  is  the  surface  height  correlation  function.  Since 
the  surface  rouqhness  is  assumed  to  be  isotropic,  the  height  cor¬ 
relation  function  R^Tx,ty)  is  a  function  only  of  the  separation 


r “  between  thg^sucface  points  x,  y  and  x',y',  (l.e.,  r  «  J{x-x')2 
+  (y-y*)  aJTv+V*  Hence  one  can  define  p(r)  «  R(t  ,t  )  to  be 
the  surface  height  correlation  function.  Here  1r  Is  th^surface 
mean  square  height  and  is  given  by  the  relation 

h2  =  <^(x,y)>. 

Thus  the  scattering  cross  section  Is  proportional  to  the  power  spec¬ 
tral  density  of  the  surface  height.  The  height  spectrum  of  the 
ocean  surface,  S(k  ,k  ),  has  been  reviewed  by  Phillips10,  who  found 
that  In  the  equilioriom  range,  the  spectrum  has  the  convenient  ap¬ 
proximate  form 

0  k  < 

S(k  k  1  =  S(k  1  =  <  (201 

BkQ  kj  <  k  <  kc, 

where  kQ  =  Jk2+kJ, 

8  Is  the  Phillips  constant,  B  -  0.006, 

kl  =  9/u  v 
kc  =  3.6  cm-1, 

g  is  the  acceleration  of  gravity  (9.8  m/sec2)  and 
is  the  wind  velocity  (m/sec). 


The  lower  cut-off  kj  represents  the  gravity  waves  with  phase  velocity 
equal  to  the  wind  velocity.  The  upper  cut-off  kc  represents  the 
capillary  waves  of  minimum  phase  velocity.  For  wave  numbers  greater 
than  i^,  the  exact  form  of  S(k0)  is  not  well  known,  and  is  very  sen¬ 
sitive  to  local  wind  speed.  Thus  at  present,  the  use  of  the  per¬ 
turbation  theory  is  limited  to  electromagnetic  wavelengths  greater 
than  about  4  cm.  Physically,  the  Phillips  spectrum  implies  that 
the  wind  speed  does  not  affect  the  shape  of  the  spectrum  in  the 
equilibrium  range.  As  the  wind  speed  increases  it  merely  drives 
the  cut-off  (k.)  to  smaller  wave  numbers.  The  upper  cut-off,  k  , 
occurs  because1viscous  dissipation  quickly  damps  out  larger  wav^ 
numbers.  Figure  9  shows  measurements  (of  Valenzuela  et  a  1 . 1 ^ )  of 
the  spectrum  of  the  sea  surface  S(k) , computed  from  measured  back- 
scattering  cross-section  at  UHF,  L,  and  C  band.  It  is  apparent 
that  the  Phillips  spectrum  is  an  excellent  approximation  to  the 
actual  spectrum  in  the  wave  number  region  where  the  perturbation 
theory  is  valid,  and  that  the  cross-sections  based  on  Equation  (19) 
and  (20)  are  in  good  agreement  with  measured  results. 


Figure  9.  Spectra  inferred  from  vertical  radar  cross-section 
data  hy .means  of  the  composite  model,  after  Valenzuela 
et  al.  .  The  line  indicates  the  Phillips  spectrum 
for  wind  speeds  of  10  and  ’0  ms. 

B.  The  Specular  Component;  Very  Rough  Surface  Scattering 

The  specular  component  a  is  calculated  from  the  "very  rough 
surface"  scattering  model.  Tnl  very  rough  class  of  surfaces  are 
those  having  a  large  scale  roughness  such  that  the  r.m.s.  rough¬ 
ness  height,  h,  is  larger  than  the  electric  wavelength  and  the  local 
surface  radii  of  curvature  are  also  greater  than  this  wavelength. 

In  the  past,  two  different  techniques  have  been  used  to  calculate 
the  fields  scattered  by  such  surfaces,  viz.  the  Physical  Optics 
approach,  and  the  Geometrical  Optics  or  stationary  phase  approach. 

If  these  are  used  correctly,  they  lead  to  the  same  result. 

The  Physical  Optics  approach,  which  employs  a  physical  optics 
integral  expression  for  the  scattered  field,  is  probably  the  most 
mathematically  thorough  and  exact  approach.  Among  others,  Semenov  , 
Kodis1J  and  Barrick1^  have  addressed  the  problem  of  scattering  from 
a  very  rough  surface  of  homogeneous  material,  using  a  vector  formu¬ 
lation  which  accounts  for  polarization.  Since  this  method  probably 
yields  the  least  physical  insight  into  the  scattering  mechanism 
and  is  computationally  time  consuming,  it  will  not  be  used  here. 

1C 

The  Geometrical  Optics  approach,  originated  by  Muhleman  , 
subdivides  the  rough  surface  into  a  grid  of  smaller  flat  sub-areas 
each  havinq  the  same  projected  area  in  the  x-y  plane.  Each  sub- 
area  reflects  the  incident  power  specularly,  and  the  direction  of 
reflection  is  determined  by  the  direction  of  the  local  normal  fi.e.. 
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the  direction  bisecting  the  angle  between  the  incidence  and  scat¬ 
tering  directions!.  Hence  the  amount  of  power  scattered  into  a 
given  direction  is  directly  proportional  to  the  number  of  sub-areas 
whose  normals  are  pointed  in  the  proper  direction,  divided  by  the 
total  number  of  sub-areas.  This  quantity  is  the  probability  density 
function  for  the  surface  normal.  An  identical  result  for  the  scat¬ 
tering  cross  section  may  also  be  obtained  from  a  slightlv  different 
point  of  view,  in  which  the  incidence  and  scattering  directions 
are  specified,  and  the  specular  reflection  points  on  the  surface 
located.  From  elementary  geometrical  optics  considerations,  the 
scattering  cross  section  of  a  single  specular  point  is  tt  times  the 
principal  radii  of  curvature  at  that  point.  For  multiple  specular 
points,  Kodis  ,  starting  with  an  exact  integral  equation  for  the 
fields  and  evaluating  it  asymptotical ly  by  the  method  of  stationary 
phase,  has  shown  that  the  average  scattered  power  is  proportional 
to  the  average  number  of  specular  points  on  the  surface  times  the 
average  radii  of  curvature  at  these  points. 

From  these  two  analyses,  Barrick^  has  presented  the  equation 
for  the  average  scattering  cross  section  per  unit  area  of  a  very 
rough  surface  of  homogeneous  material.  This  is  given  by 

°nc = "  nA<irr?i>iRr1c(v)i2 

where  n^  =  average  number  of  specular  points  per  unit  area, 

< t r*i r^ ( >=  average  absolute  value  of  product  of  principal  radii 
of  curvature  at  specular  point, 

R  ( v!  =  Fresnel  reflection  coefficient  for  incident  polarization 
ne  state  c  and  scattered  polarization  state  g.  These  quan¬ 
tities  are  given  in  Appendix  C. 

v  =  local  angle  of  incidence  at  specular  point  f half  the  angle 
between  incidence  and  scattering  directions!  with  cosv  = 

Sfl -sinPisindscos4)s+cosBicos0s,  see  Figure  10. 

After  determining,  from  the  statistical  properties  of  the  surface, 
the  quantities  n.  and  <|r.r?|>,  Barrick  obtained,  for  a  general 
surface,  ^ 

o“£  ■  »  sec4B  p(Cxsp.Cysp)  IV»)|2 

where  P(^XSD»fVSD)  is  the  joint  probability  density  function  for 
the  surfacirslopes  in  both  x,y  directions  at  the  specular  point 
with  £XSpand  r’vsD  defined  as  (9£(x,y))/3x  and  (3S(x,y)  )/3y  at  the 
specular  point*  PR  is  the  angle  between  z  axis  (mean  surface  normal) 
and  the  local  surface  normal  with  cos  fi  =  (cosP  +cos0.)/2  cosv  (see 
Figure  10).  Since  the  probability  density  function  fdr  the  surface 
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Figure  10.  Geometry  of  bistatic  cross-section. 

slop's  in  the  x  and  y  directions  is  simply  related  to  the  probability 
density  function  for  the  local  surface  normal,  one  can  show  that, 
if  the  slope  distribution  is  independent  of  azimuth  angle  4>n ,  then 

V  =  nsecB  P(0)|RnG{v)!2  (23) 

where  P(8)d,Q  is  now  the  probability  that  the  surface  normal  lies 
in  a  small  cone  of  solid  angle  d n i n  the  direction  n1  (see  Figure 
10).  The  difference  factor  sec^B  between  Equations  (22)  and  (23) 
is  due  to  the  Jacobian  of  the  transformation  between  8,  <1>  and 
C  ,  C  .  The  probability  density  function  P(8)  for  the  ocean 
surrace^depends  on  a  number  of  factors,  but  for  the  purpose  of  this 
report  it  is  assumed  to  have  the  simple  isotropic  form  (very  simi¬ 
lar  to  the  cross-wind  form  of  the  Cox-Munk16  distribution): 

P(8)  =  exp[-tan?8/t?(l+t?)j  (24) 

,  nc 

where  t  is  the  mean  square  surface  slope.  For  the  sea  surface 
which  possesses  the  Phillips  spectrum,  the  effective  mean  surface 
slope  is 


t  =  R  loq(k,,/k,)  (?R) 

Here  k1  and  k9  are  the  effective  lower  and  upper  limits  of  the  ocean 
spectrum,  respectively.  The  lower  wave  numher  limit,  in  the  Phillips 
spectrum,  is  qiven  in  Equation  (’fl).  The  upper  limit  is  determined 
hy  the  condition  that  there  he  a  smooth  transition  in  the  back-scat- 
terinq  cross  section  between  the  results  of  Equations  (?3)  and  (IQ). 
The  expression  qiven  here  is  a  convenient  empirical  result: 
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? 

(3.5X)2B 


( ?6) 


Both  shadowing  and  multiple  scattering  between  different  surface 
points  are  neglected  in  the  above  theories  for  the  very  rough  surface; 
these  limitations  are  not  expected  to  be  serious  as  long  as  the 
slopes  are  not  too  precipitous  and  near-grazing  angles  are  avoided. 

If  it  is  desired  to  account  for  the  shadowing  correction,  there 
is  a  large  literature  (see,  for  example,  reference  17)  devoted  to 
the  problem. 

C.  The  Composite  Cross-Section 

The  "composite"  model  for  the  surface  views  the  small  scale 
structure  of  the  surface  (which  produces  the  Bragq  scatter  compo¬ 
nent)  to  he  superimposed  on  the  large  scale  structure.  Thus  the 
perturbation  component  of  the  scattered  field  is  viewed  as  a  local 
scattering  phenomenon  occurrinq  on  a  tilted  surface  which  must  be 
averaged  over  the  larqe  scale  slope  distribution  of  the  surface. 

The  general  procedure  has  been  discussed  is  detail  hy  Valenzuela,  [25] 
hut  for  the  Durpose  of  this  report  it  is  more  convenient  to  consider 
expressions  for  the  bistatic  and  the  hack-scatterinq  cross-section 
separately. 

1 •  Backscatter 

For  back-scatter,  (e.=0  •  4>=n)  both  the  specular  and  the  Bragg 
components  of  scatter  will  be  significant.  When  the  averaging  of 
the  Bragg  component  over  the  slope  distribution  is  carried  out, 
the  resulting  cross-section  can  be  written  in  the  form 

o^c ( 0 i )  =  7TSec0i  lRnG( o)  | 2  P(6.)  +  16nk4cos49.  |a^|2S(?k  sin9.,0) 

(27) 

where  we  have  chosen  the  specific  polarization  states  h  (horizontal) 
and  v  (vertical)  to  represent  the  incident  polarization  (q),  and 
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that  polarization  component  ( cl  of  the  scattered  field  accepted 
hy  the  receiver  ( thus  represents  a  transmitter  with  horizontal 
polarization  and  a  receiver  with  vertical  polarization).  Clearly, 
for  hack-scatter,  the  Fresnel  coefficients  |R^o)|  ?and  l8uv^0^ 
are  equal,  while  R.  =  R  .  =  0.  The  parameter  |a'  \  includes  the 
effect  of  the  average  over  the  surface  slopes,  with  direct  polarized 
form 


| oi *  | ^  =  | at  I'*  cos^y  +  | a  sin?y  (?8a) 

1  nn  no  1  ee1  1 

where  n  =  h  or  v,  e  =  v  or  h,  n  t  e  for  and  0s=tt , 

and  cross  polarized  form 

Kel2  =  l°nel2  cos2y  +  lacnl2  si"2Y  (28b) 

where  n  =  h  or  v,  e  =  v  or  h,  n  i  e  for  9^  =  9$  and  4>s  =  -rr. 

2 

The  factor  sin  y  represents  the  effect  of  averaging  over  the  surface 
slopes  and  can  be  approximated  by  the  convenient  empirical  form 


.  P  t 
sm  y  =  Tj- 


rt^+sinfl^ 

f  t'>+sin'V):> 

i 


(?Q) 


where  t.  is  the  r.m.s.  surface  slope  given  by  Equation  (?B).  The 
coefficients  a  are  given  in  Appendix  B.  Thus  the  back-scattered 
clutter  contribution  to  the  signature  from  each  patch  is  calculated 
from  Equations  M7),  (77),  (PS)  and  (70).  Note  that  for  each  patch 
in  the  model,  the  local  scattering  anqles  9.,  9  ,  $  ,  etc.  must 
he  computed  from  the  position  of  the  radar  and  the  position  of  the 
patch  center. 


7.  Bistati'-  scatter 

Although  in  general  the  bistatic  cross-section  should  contain 
both  the  perturbation  and  specular  components  of  the  cross-section, 
it  is  found  that  in  the  forward  scatter  direction  (for  <f>  <  n/2) 
the  total  cross-section  is  dominated  by  the  specular  par!.  Thus 
in  the  ship-sea  model,  the  ship-sea  interaction  term,  which  involves 
a  forward  scatter  contribution  from  the  surface,  is  found  from  the 
specular  component  alone  (Equation  (?T)).  Again,  however,  the  local 
scattering  angles  0.,  8  and  4>  which  give  3  and  v  (and  the  appro¬ 
priate  poiarization'states)  mult  be  found  from  the  position  of  the 
radar,  the  patch  center,  and  the  specular  point  on  the  ship  structure. 
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CHAPTER  V 

DOUBLE-REFLECTION  SIGNATURE  FROM  THE  SHIP  SEA  INTERACTIONS 


A.  Numerical  Search  Technique 
for  Points  of  Reflections 


In  this  section  the  techniques  for  calculating  the  ship-sea 
interaction  term  are  considered  in  detail.  In  the  flat  patch  model 
for  the  sea  surface  described  in  Chapter  IV,  the  scattering  from 
each  patch  is  calculated  in  terms  of  the  area  of  the  patch  and  its 
bistatic  radar  cross-section  as  if  the  entire  scattered  field  em¬ 
anated  from  a  single  reference  point  in  the  patch.  Once  the  "re¬ 
flecting"  properties  of  the  patch  of  ocean  are  known,  one  may  return 
to  the  problem  of  calculating  the  double-bounce  contribution  of 
a  particular  patch.  From  any  one  of  the  sub-areas  to  the  receiver 
via  the  reflection  from  the  elliptical  cylinder,  there  must  exist 
one  reflection  point  on  the  cylinder  except  when  the  receiver  is 
in  the  shadow  reqion  of  the  cylinder.  To  find  this  reflection  point, 
the  laws  of  reflection,  as  defined  by  Equations  (?a)  and  (?b)  must 
he  satisfied.  When  these  laws  are  applied  in  a  systematic  manner, 
it  is  found1"  that  the  position  of  the  reflection  point  can  be  re¬ 
covered  from  a  sinqle  unknown  parameter  a  (defined  in  the  reference 
oivenl  which  satisfies  a  sixth  order  polynomial, 

f^a^+C^a^+C^m^+C^a^+C^oi^+C^a^  +0^=0.  f  101 

The  coefficients  C  ,  C.,  etc.  are  functions  of  the  "source"  location 
M.e.,  the  center  of  tne  patch!,  the  receiver  location,  and  the 
cylinder  parameters.  Since  the  derivation  of  this  polynomial  is 
quite  a  laborious  process,  and  it  has  been  derived  in  the  reference 
qiven  above,  it  will  not  be  reproduced. 

The  polynomial  technique  is  very  time  consuming  if  it  must 
he  carried  out  for  every  sub-area  at  each  point  of  the  trajectory. 
Hence  it  is  used  only  for  initial  acquisition  of  the  reflection 
point  for  the  first  sub-area  of  the  Nth  row  (A,...)  as  shown  in  Figure 
11.  This  reflection  point  is  then  tracked  by  irvery  fast  numeri¬ 
cal  search  technique  through  many  small  increments  (width  =  AS.) 
to  the  new  reflection  point  for  sub-area  A..,,.  This  numerical  search 
technique  is  thereafter  applied  to  find  allrthe  reflection  points 
for  the  sub-areas  in  the  same  row  (sub-area  AN-,  A,..,  etc.),  and 
the  procedure  is  then  repeated  for  each  row. 

The  polynomial  method  can  be  considered  as  an  exact  solution 
f except  for  a  numerical  truncation  error).  In  the  numerical  search 
technique,  a  new  reflection  point  is  found  from  the  previously  de¬ 
termined  point;  therefore  this  approach  can  have  an  accumulation 
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Fiqure  11.  Sub-divided  patches  of  the  illuminated  sea  surface. 


of  errors.  In  order  to  control  this  situation,  an  error  bound  is 
defined  such  that  if  it  is  exceeded  by  the  accumulated  error,  the 
numher  of  increments  is  increased  (Ail  being  reduced)  until  the  error 
bound  is  satisfied.  The  numerical  search  technique  is  briefly  dis¬ 
cussed  below.  One^needs  first  to  define  (see  Figure  3)  position 
vectors  s,  o,  and  r,_wh ere  the  vector  s  is  the  position  of  the  center 
of  sub-area,  vector  o  is  the  position  of  the  observation  point  (which 
is  fixed  for  any  instantaneous  antenna  location)  and  vector  r  is 
the  position  of  the  reflection  point  on  the  cylinder.  Thus  the 
vectors 


.  A  A 

s  *  x$x  +  y  y 


tf  - 
$  = 


V  +  V  +  V 


V  +  yry  +  Trz 


define  the  corresponding  coordinates  x^,  y  etc. 

A  parameter  v  is  introduced  to  define  the  reflection  point 

on  the  elliptic  cylinder  such  that 

x  y 

r  Jr 

cos  v  =  — -  ,  sinv=|p 

where  a,  h  are  the  semi-major  and  semi-minor  axes  of  the  ellipse. 
Thus, 

A  A  A 

"r  -  a  cos  v  x  +  h  sin  v  y  +  z  z. 

r 

Now  Equations  ("’a!  and  (?hl  can  be  combined  to  form: 

fn  •  1)  fn  x  R)  +  fn  x  Ilfn  •  Rl  *  0  (31) 

\  A  A  A 

where  n  =  n  x  +  n  y  =  b  cos  v  x  +  a  sin  v  y  a  normal  vector 

X  .7 

<t?a) 

^  A  A  ^  J, 

I  =  I  X  +  I  V  +  I  z  =  r  -  s 

a  y  Z 

*  A  A 

-  (a  cos  v  -  x  ) x  +  (b  sin  v  -  y  ly  +  z  z  an  incident 
"  r  vector 

(3?b) 

^ 

R  =  R  x  *■  R  v  +  R,z  =  0  -  r 
x  y-  z 

A  A 

=  (x  -a  cos  v)x  +  (y  -b  sin  v)y  +  (z  -z  )z  a  reflection 

0  r  vector 

(3?cl 

After  performing  the  dot  and  cross  products.  Equation  (311  can  be 
separated  into  two  equations: 


(33hl 


By  substituting  Equations  (32b)  and  (32c)  into  Equation  (33b),  one 
finds  that 


("xW^o 

'VxSV-'W'V'y' 


(34) 


which  is  a  function  of  v,  x  and  y  .  For  Equation  (33a),  one  can 
apply  the  Taylor  Series  expansion  which  says  that  the  equation  for 
the  new  reflection  point  f f .  ,  =0)  is  given  in  terms  of  the  previous 
point  (fj)  by  J 


f  J+,  .  f  J  ♦  <fj  ■  0 

3f4  3f 
where  6fjS_I6v+_ 


3f  1 

with  indicating  the  value  of  f1-  evaluated  at  the  jth  point. 

3xs  3xs 

Hence, 


3fj 

5x'  +  3^  6ys 


3f 


Sv  =  - 


fj+(3fJ./3xs)6xs  +  (3fj/3ys)6y$ 
3fj/3v 


(35) 


which  is  a  function  of  the  old  reflection  point  (v,x  ,y  )  and  the 
incremental  distances  (6x  ,6y  ).  The  new  ref lectionspo¥nt  in  the 
x-y  plane  is  now  defined  By  v^  ,  *  v.  +  6v  and  can  be  substituted 
along  with  the  new  source  location  (x  +6x  ,  y  +6y  )  into  Equation 
f 34)  to  yield  the  new  reflection  point  insthexz-coordinate.  After 
the  point  (j+1)  has  been  determined  it  can  be  substituted  into  Eq¬ 
uations  (351  and  (34)  again  to  find  the  next  new  point  (j+2).  This 
process  continues  until  it  reaches  the  center  of  the  adjacent  sub- 
area  from  which  an  additional  reflected  field  is  then  calculated. 


B.  Reflection  from  the  Elliptical  Cylinder 

Once  the  field  scattered  from  the  center  of  each  sub-area  has 
been  calculated  and  the  corresponding  reflection  point  on  the  cyl¬ 
inder  located,  the  reflected  field  incident  on  the  receiver  is  readily 
computed  using  the  geometry  illustrated  in  Figure  12.  The  scattered 
field  for  each  sub-area,  as  described  in. Chapter  IV,  is  treated 
as  the  field  incident  on  the  cylinder  (E^Qn)).  The  reflected  field 
is  given  by  Equation  (1),  K 


Er(s) 


=  E^Qu)  .  R 


r  r 
plp2 


,-jks 


(Pj+shpIJ+s) 


31 


ANTENNA 


SEA  -  PATCH 
SOURCE 


Figure  1°,  Geometry  for  field  reflected  from 
elliptical  cylinder. 

In  order  to  calcylate  the  principal  radii  of  curvature  of  the 
reflected  wavefront  p1  and  pT,  the  values  of  R..  and  R~,  the  princi 
pal  radii  of  curvature  of  the  elliptic  cylinder,  are  needed.  In 
this  case 


Using  the  above  expressions  and  Equation  (4),  one  finds  that 

12!*!  and  Vf 
Pj  s  RjCOse 

All  the  parameters  have  been  defined  before  in  Chapter  III.  The 
dyadic  reflection  coefficient  is  given  by 

■  A  ^  Am  A  A 

R  =  fe  e  -eel 

L  ii  ii  x  a  J 

and  can  he  calculated  in  the  following  manner  f see  Figure  IT). 
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Figure  13.  Reflection  coordinate  transformation  system 
on  the  elliptical  cylinder. 


The  incident  vector  I  in/igure^l3  can  be  transformed  to  the 
orthogonal  coordinate  system  n,  e^,  e^  as 

f  =  -cose’n  -  sine^osu  e?  -  sine1  sinu>  e^  (36) 

where  w  is  the  angle  between  e?  and  the  projection  of  f  onto  the 
plane  containing  e^  and  e^  sucn  that 

I  •  e?  =  -sine  cosw  . 

/V  A 

Since  it  is  also  true  that  I  •  e,,  =  -cose^,  we  have 

cose? 

cos  w= - r 

sine 

The  unit  vector  e.,  is  found  by  the  orthogonal  relationship: 


e 


1 


where  n  can  he  obtained  from  Equation  (3?a). 


el  =  e1xx  +  V 


Hence 
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wi  th 


a  sin  v 


1.2  2^2.2 
J  b  cos  v+a  sin  v 


and  e. 


b  cos  v _ 

~2  22.2 
b  cos  v+a  sin  v 


From  Figure  13  one  can  easily  transform  the  vectors  used  in  the 
dyadic  reflection  coefficient  into  the  x-y-z  coordinate  system, 

e4  =  -sinu  e^  +  cosiu  e^ 

=  -sinw  z  +  cos  to  •  e1x  x  +  cosoi  •  e^  y  (37al 

and 

e,1,  5  !  x  f4  ( 77b) 

ell  =  6  x  eA  (?7C) 

Due  to  the  principle  of  reciprocity,  the  ray  which  travels 
in  the  opposite  direction  to  that  shown  in  Figure  1°  will  produce 
exactly  the  same  field  at  the  receiver  terminals;  therefore  the 
douhly  reflected  field  calculated  for  the  ray  path  depicted  in  this 
figure  should  be  doubled  to  give  the  total  field. 


Test  Cases:  The  Elliptic  Cylinder;  The  Reflection  Points 


Before  showing  the  results  for  the  doubly  reflected  field, 
it  is  helpful  to  ensure  that  the  GTD  solution  for  the  reflected 
field  from  an  elliptic  cylinder  is  correctly  calculated.  Consider 
the  test  case  with  an  electric  current  moment  placed  close  to  an 
elliptic  cylinder,  as  shown  in  Figure  14,  with  the  reflected  field 
pattern  calculated  in  accordance  with  Equation-fl) .  The  result 
is  compared  against  that  obtained  by  Marhefka  ,  and  excellent  agree¬ 
ment  (see  Figure  15)  is  achieved.  Marhefka's  solution  which  is 
also  based  on  the  GTD  technique  has  itself  been  validated  by  a  moment 
method  solution. 


The  results  of  the  computer  program  for  obtaining  the  reflection 
points  on  the  elliptic  cylinder  have  also  been  verified.  As  an 
example,  for  the  case  of  broadside  incidence  on  the  cylinder,  as 
shown  in  Figure  IF,  the  locations  of  the  reflection  points  on  the 
cylinder  versus  the  sub-area  positions  are  illustrated  in  Figures 
1 7  and  IB.  Figure  17  shows  the  reflection  points  in  the  x-y  plane 
as  measured  by  the  parameter  v  f the  elliptic  parameter  v  is  related 
to  angle  (t>  hv  v=arctan|  (A/Bltan<f>|  1,  and  Figure  IB  shows  the  reflec¬ 
tion  points  in  the  z-plane  as  measured  by  the  height  zr.  For  the 
case  of  off-broadside  incidence  (45°),  shown  in  Figure  IB,  the  re¬ 
flection  points  versus  the  positions  of  the  sub-area  are  illustrated 
in  Figures  90  and  ?i .  The  curves  of  position  generated  for  broad¬ 
side  incidence  are  symmetric  with  respect  to  the  center  sub-area 
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Figure  14.  Configuration  for  the  reflected  field  pattern 
measurement.  The  source  is  located  in  the  x-y  plane. 


PRESENT  SOLUTION 


Figure  15.  Reflected  field  pattern  for  the  geometrical 
configuration  shown  in  Figure  14  with  A=2.X,  B=lX, 
R=4X,  4>s=45°,  RF=200X,  ZF=200\. 
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Figure  16.  Characteristic  views  of  elliptical  cylinder  and 
position  matrix  for  broadside  incidence.  The  antenna  is 
instantaneously  fixed  at  location  (x=55.2,  y=0,  z=40.5, 
and  a=3,  b»30;  arbitrary  distance  units). 
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Figure  17.  Reflection  points  (measured  in  angular  parameter  vl 
versus  the  sub-area  location  for  the  case  shown  in  Figure  16. 
The  sixth  row  in  the  scattering  matrix  is  used. 
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Figure  1°.  Characteristic  view  of  elliptical  cylinder  and 
position  matrix  for  off-broadside  incidence  of  the 
antenna.  The  antenna  is  instantaneously  fixed  at 
location  (x=3Q.n,  y=3Q.O,  z=4f).s  and  a=3,  b="*0; 
arbitrary  distance  units!. 
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Figure  90.  Reflection  points  (measured  in  angular  parameterl 
v)  versus  the  sub-area  location  for  the  case  shown  in 
Figure  19.  The  seventh  row  in  the  scattering 
matrix  is  used. 
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Fiqure  ?\.  Reflection  point  (measured  in  vertical 
heiqht  -  z  )  versus  the  sub-area  locations  for 
thersame  case  shown  by  Fiqure  ?0. 


4? 
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.is  exported.  The  curves  in  Fiqures  I7  and  1R  are  nearly  st.raiqht 
lines  for  the  reason  that  the  hroad  fare  of  the  elliptical  cvl  Intier 
is  nearly  flat.  The  abrupt  chanqe  at  the  end  of  the  curves  in  Fiq- 
ures  °0  anti  ?1  is  due  to  the  fact  that  the  reflection  points  have 
rearheti  the  narrow  end  of  the  elliptical  cylinder  where  the  surface 
curvature  varies  drastically  with  position. 

n.  Results 

After  the  above  technique  and  tests  have  been  verified,  we 
are  now  ready  to  consider  the  actual  scattered  fields  at  the  receiver 
due  to  the  double  bounce  corner  contribution  f forward  scatter  from 
transmitter  to  ocean  to  ship  to  receiver!  as  the  radar  position 
moves  alonq  a  straiqht  line.  For  each  sub-area  of  the  sea  surface 
seatterinq  grid,  the  field  incident  upon  it  is  calculated  from 
the  antenna's  radiation  functions  as  expressed  in  Appendix  A.  lhe 
scattered  field  from  this  sub-area  is  then  obtained  by  multiplyinq 
the  incident  field  bv  a  seatterinq  coefficient  which  is  a  function 
of  the  anqle  of  incidence,  anqle  of  reflection,  sea  surface  char¬ 
acteristics,  polarization  of  the  incident  field,  and  size  of  the 
sub-area  as  qi von  by  Fquations  (17!  and  (73!.  It  is  clear  that 
the  anqle  of  reflection  is  determined  bv  the  location  of  the  re¬ 
flection  point  on  the  cylinder.  This  location  is  first  calculated 
b v  the  polynomial  method  and  then  tracked  by  the  numerical  search 
technique  as  the  seatterinq  source  moves  from  one  sub-area  to  another. 
Once  the  scattered  field  from  the  sub-area  to  the  reflection  point 
on  the  cylinder  is  computed,  the  reflected  field  from  this  reflection 
point  to  the  receiver  is  readily  calculated  bv  the  method  of  GTO 
which  is  qiven  by  Tqiiation  (H.  The  final  solution  is  then  a  sum¬ 
mation  of  all  the  doubly  reflected  fields  for  the  sub-areas  that 
are  illuminated  bv  the  antenna.  Of  course,  the  field  at  the  receiver 
terminals  needs  to  be  multiplied  bv  the  antenna  pattern  function 
aqain.  A  number  of  tvnical  siqnatures  for  the  qoometrv  illustrated 
in  Tiqure  have  been  calculated,  and  are  shown  in  Fiqures  07 
throuqh  '’•> . 

In  these,  and  subsequent  siqnatures,  the  oloctromannetic  wavo- 
lenqth  has  been  taken  as  the  unit  of  lenqth;  i.e.,  the  dimensions 
of  the  ship  models,  and  the  altitude  f7F!  and  distance  alonq  the 
trajectory  (RFl  are  qiven  in  units  of  wavelenqth.  The  siqnatures 
then  represent  the  nominal  received  voltaqe  V.  in  dB  versus  the 
distance  RF  of  the  antenna  from  its  startinq  position.  Note  that 
the  receiver  voltaqe  is  calculated  only  at  finite  intervals  alonq 
the  trajectory,  the  incremental  distance  heinq  1.7  lenqth  units 
in  all  siqnatures  (except  Fiqures  66-71!.  THe  received  voltaqe 
has  been  interpolated  between  the  calculated  points  to  produce  a 
rent  1  minus  siqnature.  Note  also  that  the  vertical  scales  of  the 
plots  indicate  the  nominal  received  voltaqe  V,  in  dB  (i.e.,  V,  - 
70  loq  F^l,  where  F^  is  the  received  electric* f ield  intensity1 in 
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TOP  VIEW 


Figure  ??.  Configuration  of  elliptical  cylinder  at  sea. 
Patch  size=3x3  length  units.  Vertical  polarization 

is  employed. 
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Figure  '’t.  Voltaqe  siqnat.ure  for  the  model  of  Fiqure  r*?. 
Horizontal  trajectory  (AT=0°1  at  altitude  40. R  units. 


Ship  dimension  A=5,  B=?5  units.  Antenna  depression 
angle  ATA=60°,  Antenna  approach  anqle  SAH=0°. 

RF  is  the  distance  traveled  by  the  radar  from 
its  starting  location.  The  ship's  center  is 
located  at  RF=1?0  lenqth  units.  Cylinder 
height,  is  ?4.t6  length  units. 


volts  per  meter  calculated  by  the  equations  employed  in  this  report. 
If  one  wishes  to  calculate  the  power  delivered  (P.  , )  to  the  matched 
load  of  the  antenna  from  the  scattered-voltage  plol$  given  in  this 
report,  given  the  power  transmitted  (P.)  from  the  antenna,  the  fol¬ 
lowing  equation  is  available 


P 


del 


V-/10 

in  r 


iWeff 


ns) 


where  G.  is  the  transmitting  antenna  directive  gain,  A  „  is  the 
receiving  antenna  effective  aperture,  J  and  a  are  def^n&d  in  Ap¬ 
pendix  A,  k  is  the  free  wave  number,  and  Zo  is  the  free  space 

impedance.  Note  that  10  r  recovers  the  quantity  in  (volts/ 
meter)  . 


Also  in  computing  the  signatures,  it  has  been  assumed  that 
the  dielectric  constant  is  approximately  that  of  sea  water  at  S 
band  (e  =69-j38),  and  the  r.m.s.  slope  of  the  ocean  surface  was 
10°,  and  the  antenna  current  has  been  taken  as  J  a  -1.  In  com¬ 
puting  the  back-scattered  and  double-bounce  contributions,  the  patch 
sizes  were  taken  as  3x3  length  units  (except  in  Figures  66-71)  and 
the  patch  array  size  was  17x17  (ie.,  289  patches).  Vertical  polari¬ 
zation  of  transmitter  and  receiver  was  used,  unless  otherwise  speci¬ 
fied. 


In  Figures  23  through  25,  signatures  with  horizontal  tra¬ 
jectory  ( AT=0°)  made  with  a  15°  antenna  beam  at  a  depression  angle 
(ATA)  of  60°  are  shown,  for  the  case  of  the  elliptic  cylinder  target 
with  semi-major  axis  of  25  length  units  and  a  semi-minor  axis  of 
5  length  units.  The  trajectory  altitude  (ZF)  is  40.5  length  units, 
and  the  horizontal  range  (radar  to  target  center)  varies  from  I’O 
to  0  lenqth  units.  The  three  different  signatures  in  Figures  ?3 
through  25  correspond  to  different  antenna  approach  anqles.  which 
varied  from  broadside  incident  (SAH=n°)  to  head-on  (SAH=90°).  Al¬ 
though  the  qeneral  form  of  the  signature  is  quite  similar,  being 
strongly  influenced  by  the  assumed  Gaussian  pencil  beam  antenna 
oattern,  the  trajectory  and  the  beam  depression  angles,  the  relative 
received  voltaqe  varies  markedly  as  the  approach  angle  changes. 

Among  the  three  computer  plots,  the  received  voltage  level  is  the 
lowest  when  the  approach  angle  is  45°;  this  is  because  the  energy 
that  is  concentrated  in  the  main  beam  is  being  reflected  away  by 
the  cylinder  instead  of  being  reflected  back  to  the  receiver,  as 
happened  in  the  other  two  cases.  The  reason  that  the  head-on  sig¬ 
nature  is  lower  than  the  broad-side  signature  is  apparently  because 
of  the  higher  energy  spreading  factor  dur  to  the  larger  curvature 
in  the  bow  or  stern  region. 


The  direct  backscattered  signal  ( i . e. ,  the  "clutter  signal") 
from  each  ocean  sub-area  is  computed  using  the  formulation  presented 
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Figure  26.  Ship-sea  double-bounce  signal  plus  the  ocean  back- 
scattered  signal  for  the  geometrical  model  of  Figure  22. 

Trajectory  is  horizontal  (AT=0°)  at  altitude  (ZF)  of 
40.5  units.  Ship  dimension  is  A=5,  Ba25  units. 

Antenna  depression  angle  (ATA)  is  60  .  Antenna 
approach  angle  (SAH)  is  0°.  RF  is  the  distance 
traveled  by  the  radar  from  its  starting  location. 

The  ship's  center  is  located  at  RF=1?0  length 
units.  Cylinder's  height  is  24.36  length  units. 
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Figure  ?7.  Ship-sea  double-bounce  signal  plus  the  ocean  back 
scattered  signal  for  the  geometrical  model  of  Figure 
All  the  parameters  are  the  same  as  that  for  Figure  26 
except  the  antenna  approach  angle  SAH=45°. 
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Figure  ?8.  Ship-sea  double-bounce  signal  plus  the  ocean  back 
scattered  signal  for  the  geometrical  model  of  Figure  ??. 
All  the  parameters  are  the  same  as  that  for  Figure  ?6 
except  the  antenna  approach  angle  SAH*dO. 
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in  Chapter  IV.  A  uniformly  distributed  random  phase  for  each  sub¬ 
area  path  is  incorporated  into  the  backscattered  field.  These  ran¬ 
dom  backscattered  signals,  summed  with  the  forward  sea  scattered 
double-bounce  signals,  are  illustrated  in  Figures  through 
It  can  be  seen  from  these  figures  that,  except  for  broadside  in¬ 
cidence,  the  ship-sea  double-bounce  signals  have  mostly  disappeared 
into  the  noise-like  sea  clutter  signal,.  In  other  words,  the  sea 
clutter  signal  has  contaminated  the  signal  caused  by  the  double¬ 
bounce  ship-sea  interaction.  The  sudden  termination  of  the  signal 
near  the  end  of  the  trajectory  is  caused  by  the  fact  that  signals 
from  the  illuminated  sea  patches  he.yond  the  ship  were  not  included 
in  the  calculation. 

In  considering  the  exact  form  of  the  range  dependent  clutter 
signal  (which  would  thus  be  the  time  dependent  signal  in  a  radar 
moving  with  constant  velocity!  it  is  important  to  bear  in  mind  one 
feature  of  the  computer  program  which  implements  the  calculation. 

To  simplify  certain  geometrical  factors  in  the  calculation,  the 
grid  of  sub-areas  was  referenced  to  the  antenna  position  rather 
than  the  position  of  the  ship.  That  is,  the  entire  grid  moves  over 
the  sea  surface  along  with  the  antenna.  At  each  point  of  the  tra¬ 
jectory,  a  new  set  of  random  phase  angles  was  generated  for  each 
sub-area.  It  is  clear  that  this  procedure  will  generate  a  discon¬ 
tinuous  time  function  no  matter  how  small  the  incremental  distance 
along  the  trajectory  between  calculations.  Thus  each  point  on  the 
clutter  part  of  the  signal,  and  the  clutter  modulation  of  the  double 
bounce  return  must  he  regarded  as  a  sample  from  a  different  number 
of  an  ensemble  of  signatures  (an  attempt  to  obtain  the  doppler  spec¬ 
trum  of  a  position  of  the  signal  by  takinq  its  Fourier  transform 
would  produce  essentially  a  white  noise  spectrum  rather  than  the 
correct  spectrum!.  Note  also  that  because  the  clutter  signal  is 
the  sum  of  a  number  of  siqnals  with  random  phase,  there  is  a  finite 
probability  that  total  clutter  signal  may  be  considerably  larger 
than  the  levels  indicated  for  short  intervals  froughlv  the  recip¬ 
rocal  of  the  fade  rate!. 

It  is  a  relatively  simple  matter  to  chanqe  the  program  so  as 
to  keep  the  sub-areas  fixed  relative  to  the  ship,  and  to  assign 
a  random  phase  shift  <5  to  each  sub-area  which  would  remain  constant 
for  the  entire  trajectory.  This  would  produce  a  continuous  sig¬ 
nature,  which  should  be  a  good  approximation  to  the  correct  sig¬ 
nature  for  a  "frozen"  sea  surface,  i.e.,  for  a  surface  which  does 
not  change  significantly  during  the  time  over  which  the  signature 
is  computed.  The  doppler  spectrum,  in  turn,  should  then  be  a  good 
approximation  to  the  correct  spectrum.  A  more  difficult  problem 
would  occur  if  it  were  desired  to  take  into  account  the  effect  of 
the  actual  motion  of  the  sea  surface  on  the  time  signal  and  the 
spectrum.  In  this  case  the  <5  should  be  random  functions  of  time, 
varying  at  a  rate  compatible  With  the  surface  motion.  We  have  not 
addressed  the  problem  of  determining  a  suitable  form  (amplitude, 
spectrum,  etc.!  for  the  <S  (t). 


At  this  stage,  having  produced  a  proqram  to  compute  both  the 
forward  aru'  *>ac‘<  scatterinq  from  the  ocean  surface,  it  may  be  ap¬ 
propriate  to  ask  the  question:  What  is  the  appropriate  size  for 
the  suh-area  or  patch  in  the  model?  The  ocean  surface  model  is 
a  statistical  one,  with  a  uniformly  distributed  random-phase  as- 
siqned  to  each  patch.  It  is  clear  that  the  r.m. s.  scattered  field  av- 
eraqed  over  a  number  of  trajectories  should  be  independent  of  the 
size  of  each  suh-area,  as  long  as  the  total  illuminated  area  is 
kept  constant.  Even  though  Equation  (171  indicates  that  the  scat¬ 
tered  field  from  a  single  suh-area  does  depend  on  the  physical  area 
-A  ,  one  should  realize  that  when  A  is  reduced,  the  total  number 
of  sub-areas  will  consequently  be  increased.  However,  the  size 
of  each  sub-area  can  not  be  chosen  arbitrarily.  The  size  can  not 
be  too  small  otherwise  computer  time  will  become  a  problem,  and 
from  a  physical  point  of  view  the  model  will  be  incorrect  because 
it  will  assign  uncorrelated  returns  (because  of  the  random  phase 
term)  to  scattering  from  arbitrary  close  surface  points.  Nor  can 
the  patch  size  be  chosen  too  large,  otherwise  one  may  miss  some 
of  the  dominant  specular  reflection  points.  In  Figure  29  ,  a  number 
of  ship-sea  double-bounce  signatures  have  been  calculated  as  the 
sub-area  width  (pw)  was  varied  from  one  length  unit  to  ten  length 
units.  The  signatures  are  calculated  for  the  same  geometry  as  that 
for  Figures  1?  and  13  except  that  the  antenna  beam  depression  angle 
is  chosen  to  be  .TO0  instead  of  60°.  It  can  be  seen  that,  for  the 
frequency  ranqe  and  target  size  being  used,  the  average  signatures 
are  nearly  invariant  to  the  sub-area  size.  This  empirical  finding 
has  permitted  us  to  use  fairly  large  patch  sizes  of  six  to  ten  length 
units,  and  thus  achieve  reasonable  program  run  times. 
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Figure  29.  Ship-sea  double  bounce  signatures.  RF  is  the 
distance  traveled  by  the  radar  from  its  starting 
position.  Patch  width  PW  is  in  length  units. 
Antenna  beam  depression  angle  is  30.  Antenna 
approach  angle  SAH=0°.  Ship's  center  is 
located  at  RF=120  length  units. 
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CHAPTER  VI 

SHIP  SCATTERED  SIGNATURES 


Having  illustrated  the  calculation  of  the  hack-scattered  (clut¬ 
ter)  signal  from  the  sea  and  the  ship-sea  Interaction  term,  the 
scattering  from  the  ship  itself  must  now  he  considered.  In  this 
chapter  the  simplified  ship  model  Is  represented  by  two  superimposed 
finite  elliptic  cylinders,  and  the  computationally  efficient  GTD 
technique  is  utilized  to  compute  the  hackscattered  signals.  There 
are  two  principal  scattering  mechanisms  that  contribute  to  the  sig¬ 
nature  from  the  ship.  These  are  the  signatures  caused  hy  the  corner 
reflection  hetween  the  two  cylinders,  and  the  signatures  due  to 
diffracted-then-ref looted  (and  ref lected-t.hen-dif fract.ed)  fields 
hetween  the  edge  of  one  cyl inder  and  the  flat  "deck"  region.  There 
is  also  a  third  hackscattering  mechanism,  which  is  the  direct  dif¬ 
fraction  from  the  curved  wedges  of  either  of  the  two  cylinders. 
However,  since  the  back  diffracted  direction  is  not  anywhere  close 
to  the  shadow  boundaries,  as  illustrated  in  Fiqure  30,  the  direct, 
diffracted  energy  is  too  small  to  he  considered  here. 
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Fiqure  30.  Direct  wedge  diffractions  showing 
the  shadow  boundaries. 
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The  two  major  ship  scattering  mechanisms  are  described  helow: 


A.  Corner  reflected  siqnature 

The  corner  reflected  siqnature  as  illustrated  in  Fiqure  31a 
is  the  siqnature  bounced  between  the  top  flat  surface  of  the  bottom 
cylinder  and  the  curved  surface  of  the  top  cylinder.  It  is  simple 
to  show,  usinq  imaqe  theory,  that  if  the  surfaces  are  perpendicular 
to  each  other,  then  there  exists  only  one  reflection  point,  located 
at  the  corner  between  the  two  cylinders  as  shown  in  Figure  b. 

Thus,  one  is  able  to  move  the  source  point  to  its  image  location 
as  shown  in  Figure  .He,  and  the  corner  reflected  field  can  then 
be  computed  as  a  single  reflection  from  an  elliptic  cylinder.  On 
this  cylinder,  the  reflection  point  in  the  z-direction  is  always 
fixed  at.  the  height  of  the  corner;  but  the  reflection  point  in  the 
x-y  plane,  designated  by  the  angle  as  shown  in  Figure  3?,  remains 
to  be  determined.  The  elliptic  parameter  (v)  is  related  to  <j>  by 
v  =  arc tan [(Aj/B?) tan^J . 

Trun  the  laws  of  reflection  as  described  by  Equation  ('’l  and 
with  the  aid  of  the  numerical  search  technique  presented  in  Chapter 
V,  the  reflection  point  or  the  parameter  v  can  be  accurately  specified. 
Knowinq  the  location  of  the  reflection  point.  Equation  (11  can  then 
be  employed  to  calculate  the  corner  reflected  field  received  by 
the  antenna.  This  corner  reflected  field  for  the  geometry  shown 
in  Figure  33  is  plotted  aqainst  the  distance  traveled  by  the  antenna 
in  Fiqure  34.  In  this  figure,  the  three  different  curves  are  dis¬ 
tinguished  hv  the  three  different  antenna  approach  angles  (SAHl. 

The  situations  encountered  here  are  similar  to  those  observed  with 
the  ship-sea  double-bounce  signatures,  wherg  the  received  siqnal 
is  the  lowest  when  the  approach  angle  is  as1  t  and  the  head-on  re¬ 
flection  is  lower  than  that  for  broadside  incidence,  for  the  same 
reason  as  qiven  previously.  However,  one  additional  feature  that 
can  be  more  clearly  noticed  in  Figure  34  is  that  the  peak  of  the 
head-on  signal  occurs  much  earlier  than  the  other  two  due  to  the 
fact  that  the  how  of  the  ship  is  much  closer  to  the  relative  antenna 
position  than  the  broadside  face  of  the  ship.  Since  the  ship  model 
is  assumed  to  be  a  perfect  conductor,  the  above  corner  reflected 
signatures  are  expected  to  have  the  highest  voltage  level  at  the  re¬ 
ceiver  terminal  among  all  the  signatures  that  are  investigated  in 
this  report. 

B.  Pitt  >\v.ted-ref  lected  signatures 

The  diffracted-ref lected  signature  represents  the  field  that 
is  diffracted  by  the  top  wedge  of  the  bottom  cylinder  and  then  re¬ 
flected  off  the  curved  surface  of  the  top  cylinder.  By  the  prin¬ 
ciple  of  reciprocity,  the  diffracted-ref lected  signature  depicted 
in  Figure  35a  will  produce  the  same  voltage  at  the  receiver  terminal 
as  the  ref lected-diffracted  signature  shown  in  Figure  35b.  Hence 
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Figure  31.  Corner  reflection  configurations. 


Figure  32.  Top  view  of  the  ship's  corner  reflection. 
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Fiqure  34.  Corner  reflected  signatures  between  the  two  cylinders 
RF  is  the  distance  traveled  by  the  radar  from  its  starting 
locations.  The  vertical  axis  is  the  received  voltage 
level  measured  in  decibels.  Ship's  center  is  located 
at  RF=120  length  units.  Ship  dimensions  are:  Al=5, 

Bl=25,  A2=4,  B2=20  length  units.  Two  cylinders' 
heights  are  ZCY=24.36,  ZCY2=33  length  units. 
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(a)  (b) 

Figure  35.  Ship  scattering  mechanisms  for  (a)  diffracted- 
reflected  signature  (b)  ref lected-dif fracted  signature. 


the  received  field  due  to  these  two  scattering  mechanisms  is  just 
twice  the  field  predicted  by  either  one.  In  order  to  calculate 
this  term,  it  is  essential  that  one  accurately  locate  the  diffrac¬ 
tion  point  (4).)  alonq  the  bottom  cylinder's  wedge,  and  the  reflection 
point  p , z  )  on  the  curved  surface  of  the  top  cylinder  as  illustrated 
by  Fiqure  lE.  A  numerical  search  technique  is  employed  again  to 
find  these  points.  This  time  not  only  do  the  laws  of  reflection 
need  to  be  satisfied,  but  it  is  also  required  to  satisfy  the  law 
of  diffraction,  as  given  by 


e 

fel 


cos  ft 


( 39a  1 


where  e  and  p  have  been  defined  before  in  section  (111-2)  and  (III- 
3)  with  0 


e  =  -Alsin  v^x  +  Ricos  v^y  . 

-1  A1 

where  v^  =  tan  {^y  tan<j> ^ ) . 

The  incident  vector  is  given  by 

1  =  rxx  +  y  +  izz 

=  (Alcos  vt-xo)x  +  (Blsin  vryo)y  -r  (zrf-zo)z 

where  (x  ,  y  ,  z  )  is  the  antenna  location  and  z.  is  the  height 
of  the  hottom  cylinder.  The  diffracted  vector  is  given  by 
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D  =  D^x  +  Dy9  +  Dzz 

=  (A2cos  v^-Alcos  v^)X  +  (82sin  v?-8Isin  Vj)£+(zr-zd)2 
where  v?  =  tan-1  (  tan^), 

Al,  A?,  Bl,  B2  (semi -mi nor  and  semi -major  axes  of  the  lower  and 
upper  cylinders)  and  zr  is  the  height  of  the  reflection  point  from 
the  sea  surface. 

The  reflected  ray  unit  vector  direction  is  given  by 
R  =  Rx3  +  RyV  +  Rz2 

=  (xQ-A?cos  v?)x  +  (y0-B?sin  v?>y  +  (zQ-zr)z 

From  Chapter  III,  the  laws  of  reflection  can  be  combined  to  form 
fn  •  D)(n  x  R)  +  (n  x  D)  (n  .  if)  =  0  (3°b) 

where  n  =  nxx  +  n^y 

=  B2  cos  v^x  +  A2sin  v^y  . 
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Equation  (39a)  can  be  rearranged  to  give 


f/Wzr>  *  ffj 


e  - 


D 

itfl 


(40) 


and  Equation  (39b)  can  be  separated  into  two  equations  (as  in  Eq¬ 
uations  (33a)  and  (33b))  in  the  form 

f?(vi-v2>  ■  KV\V<\Rx-\V*KV\V<%VnvDx,“° 


)  =  (n  0  +n  D  )R  +(n  R  +n  R  )0  =0 
3  1  ’  ?’  r  xx  y  y  z  x  x  y  y  z 


(41) 

(4?) 


The  method  of  Taylor  Series  expansion,  as  described  in  Section  V- 
A,  can  again  be  applied  here  to  Equations  (41),  (4?)  and  (43)  to 
form  the  search  equation  given  by 
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The  new  diffraction  and  reflection  points  are  then  given  by 
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One  sees  that  Equation  (43)  implies  that,  as  the  antenna  moves 
by  an  incremental  distance  of  (<Sx  ,Sy  ,<Sz  ),  the  diffraction  point 
moves  by  an  incremental  amount  5v?  ano  the  reflection  point  by  in¬ 
crements  of  *sv~  and  6z  .  The  rignt  hand  side  of  Equation  (43)  is 
not  only  a  function  ofrthe  change  in  antenna  position  but  is  also 
a  function  of  the  previous  locations  of  (40,  (4>.>,z  )  and  (x  ,y  ,z  ). 
Therefore,  there  still  remains  the  task  of binding  the  initial  scat¬ 
tering  points  (v.l  and  (v?,z  )  so  that  the  search  equation  can  beqin 
to  generate  the  new  points.  Consequently,  one  must  start  with  the 
antenna  at  such  a  position  that  the  scattering  points  can  be  easilv 
found.  For  example,  the  antenna  can  be  initially  located  in  the 
x-z  plane  as  illustrated  in  Figure  37,  and  one  can  immediately  see 
that  v,  =  v,  =  0°.  The  reflection  point  zr  can  also  be  easilv  cal¬ 
culated  and  is  given  by 


where  all  the  parameters  are  defined  in  Figure  37. 


SIDE  VIEW 


Figure  37.  Configuration  of  the  initial  scattering  points 
for  the  search  equation. 

At  any  location  along  the  trajectory,  once  the  scattering  points 
are  found  by  the  numerical  search  technique,  the  diffracted  and 
reflected  fields  can  be  obtained  bv  Equations  (1?)  and  (1),  respec¬ 
tively.  Figure  38  shows  the  diffracted-ref lected  fields  plotted 
against  the  incremental  distance  traveled  by  the  antenna  for  three 
different  antenna  approach  angles.  The  geometrical  configuration 
used  for  these  plots  is  illustrated  in  Figure  33.  The  general  shapes 
of  these  plots  are  very  similar  to  those  obtained  for  the  corner 
reflected  field  given  in  the  previous  section.  This  is  expected, 
since  the  antenna  provides  a  concentrated  beam  pattern  and  the  size 
of  the  two  elliptic  cylinders  are  very  much  the  same,  thereby  causing 
the  scattering  points  of  the  two  scattering  mechanisms  to  be  located 
close  to  each  other.  The  amplitude  of  the  diffracted-ref lected 
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Figure  38.  Diffracted-ref lected  fields  between  two  elliptic 
cylinders.  RF  is  the  distance  traveled  by  the  radar  from 
its  starting  location.  Ship's  center  is  located  at 
RF=120  length  units.  Antenna-beam  depression 
angle  ATA=60°.  Ship  dimensions  are:  AT  =5 , 

Bl=?5,  A2=4,  B2=20  length  units.  Two 
cylinder's  heights  are  ZCY=2d.3f>, 

ACY?=33  length  units. 
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field  is  lower  than  that  of  the  corner  reflected  field  because  in 
qeneral  a  diffracted  field  is  lower  than  a  pure  reflected  field. 

Figure  .10  shows  another  diffracted-ref lected  mechanism  that 
should  be  included  in  the  calculation  for  the  total  field.  This 
scattering  mechanism,  depending  on  the  geometry  of  the  two  elliptic 
cylinders,  has  the  possibility  of  producing  fields  of  about  the 
same  magnitude  as  those  shown  for  Figure  IF.  The  reflected  field 
from  the  flat  top  (deck)  of  the  bottom  cylinder  can  be  easily  cal¬ 
culated  by  using  imaqe  theory,  and  the  diffracted  field  off  the 
top  cylinder's  curved  wedge  can  again  be  calculated  by  Equation 
(1?).  It  can  be  shown  that,  for  the  same  antenna  location,  the 
diffraction  point  on  the  top  cylinder's  curved  wedge  has  the  same 
angular  position  -  4>  as  that  obtained  for  the  corner  reflected  field 
-^shown  in  Figure  32.  For  the  scattering  geometry  shown  in  Figure 
39,  the  field  versus  the  distance  increments  traveled  by  the  antenna 
is  shown  in  Figure  40.  Since  for  the  particular  dimensions  used 
here  the  top  cylinder  is  too  high  for  the  narrow  antenna  beam  to 
cover  the  diffraction  and  reflection  points  at  the  same  instant, 
the  field  has  relatively  low  amplitude.  When  the  height  of  the 
top  cylinder  is  reduced  from  33  length  units  to  26  length  units, 
the  field  is  calculated  again  and  is  shown  in  Figure  41.  This  time 
the  field  has  a  much  higher  level  and  is  comparable  to  that  of  Fig¬ 
ure  38. 


ANTENNA 


Figure  39.  Diffracted-ref lected  field  configuration. 
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CHAPTER  VII 

COMBINED  SOLUTIONS  WITH  A  CONICAL  SHIP  HULL 
AND  DIFFERENT  ANTENNA  PATTERNS 


Thus  far  the  basic  scatterinq  mechanisms  of  the  ship  and  sea 
have  been  analyzed  individually,  and  it  has  been  found  that  there 
are  three  dominant  scattered  signatures  from  the  ship-sea  model. 

It  is  worthwhile,  at  this  point,  to  briefly  summarize  them  as  fol¬ 
lows: 

1.  Backsrattered  siqnatures  from  the  sea  surface:  Both  the 
gravity  waves  and  capillary  waves  may  contribute  signifi¬ 
cantly  to  the  hackscattered  signature. 

?.  Double-bounce  signature  between  the  ship  and  sea:  The 
forward  scattered  (bistatic)  signature  from  the  sea  sur¬ 
face  is  predominantly  governed  by  the  scatterinq  from 
gravity  waves.  The  reflected  field  from  the  ship's  hull 
(elliptic  cylinder)  is  calculated  from  Equation  (1). 

3.  Ship  scattered  signatures:  There  are  two  major  contribu¬ 
tions  from  the  two-elliptic  cylinders.  These  are  the 
corner  reflected  and  diffracted-ref lected  signatures. 

By  summinq  the  above  siqnatures  or  simply  by  adding  together 
the  siqnatures  (with  proper  phase  addition)  in  Figures  ?F  through 
with  those  in  Fiqures  34,  3R  and  40,  one  obtains  the  total  re¬ 
ceived  fields  produced  by  the  ship-sea  model,  which  are  illustrated 
in  Figures  4?,  41  and  44.  The  geometrical  conf iquration  of  the 
ship  and  ♦'he  antenna  trajectory  is  shown  in  Fiqure  3?,  For  broadside 
inciden,  a  (SAH=D°),  the  signature  as  shown  in  Figure  4?  has  a  very 
strong  ship-scattered  term  which  is  about  30  dB  above  the  average 
sea  ciutter  signal.  For  the  same  trajectory  incident  anale,  the 
ship-sea  double-bounce  siqnal  (which  occurs  at  an  earlier  time  than 
the  ship  scattered  siqnal)  has  a  peak  level  of  about  10  dB  above 
the  sea  clutter.  When  the  antenna  has  an  approach  angle  of  4F° 
as  shown  in  Figure  43,  the  ship  scattered  and  ship-sea  double-bounce 
signals  have  both  disappeared  into  the  sea  clutter  signal.  For 
the  head-on  (SAH=90°)  signature  as  shown  in  Figure  44,  the  ship 
scattered  siqnal  has  again  become  very  strong,  but  the  ship-sea 
double  bounce  signal  is  barely  detectable,  being  only  a  few  decibels 
above  the  sea  clutter.  Since  all  the  above  behavior  has  been  in¬ 
terpreted  previously  in  section  (V-D)  and  (VI-A),  no  further  com¬ 
ment  is  required  here. 

In  this  chapter  we  consider  a  further  elaboration  of  the  ship 
model  in  which  the  hull  structure  is  modified  from  an  elliptic  cy¬ 
linder  to  the  physically  more  realistic  shape  of  an  elliptic  cone. 
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Figure  43.  Ship-sea  scattered  signature  for  the  same  geometrical 
configuration  as  that  for  Figure  4?  except  that  SAH=45°. 
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Figure  44.  Ship-sea  scattered  signature  for  the  same  geometrical 
configuration  as  that  for  Figure  4?  except  that  SAH=90°. 


In  this  form,  the  model  is  versatile  enough  to  simulate  the  major 
scattering  contributors  of  wide  variety  of  ships.  In  addition, 
we  wish  to  use  a  different  antenna  radiating  characteristic  so  that 
the  effect  of  antenna  pattern  on  the  signatures  can  also  be  studied. 
The  geometrical  configuration  of  the  antenna  trajectory  and  the 
new  ship  model  resulting  from  the  modification  of  the  hull  struc¬ 
ture  is  shown  in  Figure  45.  It  is  clear  that  the  only  scattered 
signal  to  be  re-calculated  due  to  this  modification  is  the  double 
bounce  ship-sea  interaction  term.  Consequently  the  numerical  search 
technique  for  finding  the  reflection  points  must  be  modified  and 
the  GTD  calculation  for  the  reflected  field  from  an  elliptic  cone 
must  be  obtained. 

It  is  interesting  to  see  first  how  the  numerical  search  technique 
for  an  elliptic  cone  differs  from  that  for  an  elliptic  cylinder. 

For  a  general  elliptic  cone  as  shown  in-Figure  46,  the  locus  of 
its  surface  is  described  by  the  vector1^ 

A  A  A 

r  =  (z+k)tan0.|Cos  v  x  +  (z+kltanfl,,  sin  v  y  +  z  z  (44) 

where  v  is  related  to  the  <p  angle  in  the  x-y  plane  measured  counter¬ 
clockwise  from  the  x-axis;  6,,  0-  and  k  are  illustrated  in  Figure 
46.  (Note  that  e.  and  d~  are  different  from  those  used  in  previous 
chapters  and  k  is  not  the  electromagnetic  wave  number). 

The  two  principal  tangents  on  the  surface  are,  with  r  defined 
in  Equation  (44), 


rz  =  -— ■  =  tanQjCOs  v  x  +  tanS^,  sin  v  y  +  z  (45a) 

A 

ry  =  -|£=  -(z+kHanejeosv  x  +  (z+k)tane2  sin  v  y  (45b) 

The  normal  to  the  surface  is 

n  =  rv  x  rz  =  nxx  +  fiyy  +  nzz 

*  tane^cosv  x  +  tane^sinv  y  -  tane1tane2z  (46) 

The  source  location  is  on  the  surface  of  the  sea  and  can  be 
denoted  by  (x  ,y  ).  The  observation  point  at  the  receiver  is  located 
at  (x0»y0»z0)*  Therefore,  the  incident  vector  can  be  given  by 

^  A  A  A 

i  *  v +  y +  v 

=  [(Zp+kltanfl^os  v-xs] x+ PzR+kltan9?sin  v  -  yg]y  +  zR  z, 

(47) 
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Figure  46.  Elliptic  cone  configuration, 
and  the  reflection  vector  is  given  by 
R  =  Ryx  +  R^y  +  Rzz 

3  [xo-(Vk)tan0lcos  v]x+[y0-(zR+k)tan92s1n  v]y+(z0-zR)z 

(48) 

where  (z  v)  are  the  height  and  elliptic  parameter,  respectively, 
definingKthe  reflection  point  on  the  elliptic  cone.  Again  it  is 
essential  to  satisfy  the  laws  of  reflection  as  given  by 

(n  •  I)(n  x  R)  +  (n  x  I)(n  •  R)=  0  (4P) 

where  n  •  f  •  Vx  ♦  hyIy  * 

I*'/  '"yVzV1*  *  '"zVx'x'rt ^ 

n  x  R  =  (nyRz-nzRy)x  +  (nzRx-nxRz)y  +  ^Ry'W* 
n  •  R  =  nxRx  +  nyRy  +  nzRz 
Equation  (49)  can  thus  be  written  as 

ff1 +  (f?)y  +  (f3)z  =  0 
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with 


fl  *  (nxRxt"yRy+n2RZ)("y,2-nzI,)t("xIxt"y,y*n2Iz>(n»Rz-"zRy,”° 
f2  ■  <nxR«t"yVnzRz»nz,x-"xIz)t("x'xtny,ytnzIzH,,zRx-nxRz>>0 


f3  =  lnA"/J*”A",.,y-V.w".,.,Vy*nA)l,.V'Aw 


(SO) 

After  substitutinq  Equations  (46),  (47)  and  (481  into  (501,  one 
finds  that  f.,  f?  and  f3  are  all  functions  of  (v,Zg),  (x.,ys)  and 
(x  ,y  ,z  1  with  v  and  z«  as  the  two  unknown  variables,  since  there 
ar£  tRreR  quations  withKonly  two  unknowns,  one  need  only  be  con¬ 
cerned  with  two  out  of  the  three  equations.  However,  in  order  to 
ensure  that  all  the  aspects  of  the  reflection  process  have  been 
considered,  it  is  helpful  to  combine  two  of  the  three  equations 
(since  these  are  linearly  dependent  functions).  Thus  it  is  arbi¬ 
trarily  chosen  to  have 


(51) 


Once  again,  as  shown  in  sections  (V-A)  and  (VI-8),  the  Taylor  Series 
expansion  can  be  applied  to  Equation  (51)  with  the  result  given 
by 


BE  =  -f3  - 


af3  3f3  3f3  3f3  3f, 


Equation  (5?)  enables  one  to  calculate  the  shift  of  the  reflection 
point  (Sv,6zn)  resulting  from  the  incremental  distance  change  of 
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the  source  point  («x  ,6y  )  or  the  observation  point  (<5x  ,6y  ,6z  ). 
Once  again,  before  employing  Equation  (52),  one  needs  to  find  tne 
initial  reflection  point  (v,z„)  in  order  to  start  the  search.  In 
Chapter  V,  the  polynomial  technique  was  used  to  find  this  initial 
point  for  the  elliptic  cylinder.  In  the  case  of  an  elliptic  cone, 
it  is  found  that  the  polynomial  technique  is  much  too  complicated 
since  it  involves  solving  two  equations  with  two  unknowns.  However, 
if  one  places  both  observation  and  source  points  on  the  x-axis  as 
shown  in  Figure  47,  it  can  be  Immediately  determined  that  angle 
v  is  equal  to  zero  and  that  the  location  zR  is  given  by  the  equation 

zR  =  (x'-blsinftjcosfl.,. 


A 

z 


Fiqure  47.  Initial  location  of  (v,zR)  on  the  elliptic  cone  for 
the  numerical  search  technique. 

Once  the  reflection  point  on  the  elliptic  cone  has  been  initialized, 
it  can  be  tracked  along  the  surface  to  that  reflection  point  which 
corresponds  to  the  actual  locations  of  the  source  and  observation 
points,  and  thereby  the  reflected  field  can  also  be  computed.  The 
same  equations  that  were  used  to  calculate  the  reflected  field  from 
an  elliptic  cylinder  can  again  be  used  here  for  the  elliptic  cone. 
The  only  difference  is  that  one  of  the  principal  radii  of  curvature 
of  the  elliptic  cone  needs  to  be  re-calculated.  The  principal  radii 
of  curvature  are  the  surface  radii  which  lie  in  the  principal  planes 
that  contain  the  vector  normal  to  the  surface.  In  the  case  of  an 
elliptic  cone,  one  of  the  principal  radii  is  infinite,  and  the  other 
is  given  by™ 
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E+?F(feH/3v)  +  G(  azR/3v)2 

RG  - - - - - - 5—  (53) 

e+2f(azR/3v)  +  g(azR/3v) 

where 


E  = 


Ok 


F  = 


r_ 


G  = 


(54) 


here  r  and  r  are  the  two  principal  tangents  on  the  surface  and 
are  given  by  Equation  (45), 


e 


f  = 


9  = 


n 

n 

A 

n 


(55) 


with 


vv 


zv 


zz 


a2r 

3v2 


3^r 

3z3v 

3z? 


where  r  is  given  by  Equation  (44). 

If  38.  is  the  incremental  distance  on  the  surface  in  the  x-y 
plane,  one  can  show  that 


3zp 

3v 


3Zp 

38 


38 

3v 


7  ? 

sinv  cosvftan  P|-tan  p^) 


I  7  7  7  7  7  7  7  7  7  7 

>/tan  BjSin  v+tan  0j,cos  v+tan  e^tan  0?(tan  e^sin  v+tan  fl^cos  v+2) 


•7( 


z+k)2tan20^sin2v+(z+k)2tan20;,cos2v 


(56) 
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After  substituting  the  results  from  Equations  (54),  (55)  and  (5fl 
into  Equation  (53),  one  is  able  to  obtain  the  radius  of  curvature 
and,  in  turn,  the  reflected  field  from  the  elliptic  cone  using  Eq¬ 
uation  (1).  The  validation  of  this  calculated  field  is  done  by 
an  experimental  comparison  as  shown  in  Appendix  D.  A  number  of 
signatures  for  the  geometry  illustrated  in  Figure  45  have  been  com¬ 
puted  and  are  shown  in  Figures  4Q,  50  and  51.  All  the  dominant 
scattering  mechanisms  are  included.  It  appears  that  no  significant 
difference  exists  between  the  signatures  for  the  elliptic  cylinder 
and  those  for  the  elliptic  cone. 

One  additional  feature,  besides  the  effects  of  the  modified 
hull  structure,  that  can  be  observed  in  these  plots  is  that  the 
ocean  backscattered  signature  has  been  programmed  beyond  the  ship's 
center.  In  other  words,  in  the  previous  plots,  all  the  signatures 
were  calculate  for  an  antenna  position  that  starts  at  120  length 
units  from  the  ship's  center  and  ends  at  the  center  of  the  ship. 
(However  the  sea  scattered  contribution  was  calculated  only  from 
the  sub-areas  in  front  of  the  ship.)  Now,  for  the  sake  of  complete¬ 
ness,  the  antenna  starts  to  collect  signature  when  it  is  100  units 
from  the  ship's  center  and  finishes  it  at  20  units  beyond  the  ship's 
center,  with  backscattered  contribution  from  the  sub-areas  that 
are  behind  the  ship  (see  Figure  48). 


ANTENNA  ANTENNA 


Figure  48.  Configuration  for  sea  scattered  signature 
collected  beyond  the  ship. 
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Figure  49.  Total  scattered  signature  for  the  ship  with  a  conical 
hull.  Signature  is  collected  beyond  the  ship.  Geometrical 
configuration  is  shown  in  Figure  45  with  e,=4.7°  and 
e?=22.3  .  Antenna  approach  anale  SAH=0°.  1  Antenna 
“beam  depression  angle  ATA=60  .  Ship's  center  is 
located  at  RF=100  length  units  with  RF  being  the 
distance  traveled  by  the  antenna  from  its  starting 
location.  Ship's  dimensions  are:  Al=5,  Bl=25,  A2=4, 

B2=20,  ZCY=24.36,  ZCY2=33  length  units. 
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The  antenna  pattern  can  also  play  on  important  role  in  the 
form  of  the  signature.  Previously,  the  antenna  pattern  used  was 
a  15-degree  pencil  beam,  thus  confininq  all  the  significant  scat¬ 
tered  power  to  be  returned  from  only  the  small  area  illuminated 
by  the  narrow  beam.  It  is  interesting  to  see  how  the  signatures 
chanqe  when  a  qreatly  broadened  antenna  beam  is  used.  The  new  antenna 
radiating  characteristic,  as  shown  in  Figure  5?,  has  also  a  15-degree 
Gaussian  beam  but  it  is  transmitted  conically  around  the  antenna's 
trajectory.  Thus  looking  down  onto  the  ocean  surface,  the  illu¬ 
minated  area  has  been  greatly  increased  transverse  to  the  trajectory 
direction.  The  analysis  for  the  antenna  radiation  pattern  must 
now  be  correspondingly  changed. 


A 


Fiqure  5?.  Conical  antenna  pattern. 

A 

„  If  one  denotes  I  to  be  the  arbitrary  incident  unit  vector  and 

T  to  be  the  antenna  trajectory  direction,  then  the  pattern  angle 
0  from  the  trajectory  direction  is 

0  =  arccos(I  •  T)  (57) 

If  ATA  is  the  pattern  cone  half  angle,  then  the  far  field  pattern 
is  assumed  to  have  the  form 

[ka  cos(0-ATA)j? 

E(9)  =  V2  e"jkr  e_  (58) 

where  "a"  determines  the  beamwidth  of  the  pattern,  and  J  a^=i  is 
assumed  in  the  computation. 
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When  the  elliptic-cone  ship  model  is  traversed  with  the  new  hollow- 
cone  antenna  pattern,  one  obtains  the  total  ship-sea  scattered  sig¬ 
nature  shown  in  Fiqures  53,  54  and  55  for  a  cone  half  angle  (ATA) 
of  60°.  Signatures  for  the  cone  half  angle  of  30°  are  shown  in 
Figures  56,  57  and  58.  When  comparing  these  signatures  with  those 
previously  obtained  for  the  antenna  with  a  pencil  beam,  one  signi¬ 
ficant  difference  can  be  observed.  That  is,  when  the  pencil  beam 
is  used,  the  ship  can  not  be  detected  for  approach  angles  in  the 
vicinity  of  45°;  however,  when  the  hollow-cone  pattern  is  used, 
the  ship  contributes  an  observable  backscattered  signal  no  matter 
what  may  be  the  antenna  approach  angle.  This  is  true  for  both  the 
30°  and  60°  cone  half  angles.  Apparently,  the  beam  is  broad  enough 
to  pick  up  the  major  specular  points  on  the  ship  even  though  it 
is  not  normally  incident.  Another  noticeable  difference  is  that 
the  sea  clutter  signal  for  the  conical  antenna  pattern  is  higher 
than  that  for  the  pencil  beam  pattern.  This  is  true  because  more 
of  the  sea  surface  is  being  illuminated  by  the  broadened  antenna 
pattern.  An  additional  anticipated  feature  observed  in  these  plots 
is  that  when  the  cone  half  angle  is  3f)°,  the  ship  and  ship-sea  scat¬ 
tered  signal  occur  much  earlier  than  for  the  50°  cone  half  angle. 

This  is  because  the  main  antenna  beam  reaches  the  ship  much  earlier 
when  the  anqle  ATA  is  30°.  The  sea  back-scattered  signal  is  sub¬ 
stantially  lower  when  angle  ATA  is  30°,  not  only  because  the  ef¬ 
fective  range  has  been  increased  but  also  because,  at  this  anqle, 
the  backscatterinq  cross  section  of  the  ocean  decreases  as  the  angle 
of  incidence  (from  vertical!  increases,  particularly  for  horizontal 
polarization.  The  signatures  illustrated  in  Figures  53  through 
^8  are  for  vertical  polarization;  keeping  the  same  geometrical  con¬ 
figurations,  the  signatures  for  horizontal  polarization  are  illus¬ 
trated  in  Fiqures  5Q  through  54. 

Before  concluding  the  analysis  on  the  simplified  ship-sea  model, 
it  is  appropriate  to  repeat  the  calculation  with  a  model  that  closely 
resembles  a  real  ship.  The  ship  U.  S.  S.  Missouri  is  chosen,  and 
a  model  of  it  is  shown  in  Figure  65a.  The  simplified  model  with 
corresponding  dimensions  is  illustrated  in  Figure  65b.  The  calcu¬ 
lated  signatures  for  this  model  with  the  same  trajectory  and  geo¬ 
metrical  configuration  as  that  illustrated  in  Figure  45  are  shown 
in  Figures  66  through  71.  The  trajectory  altitude  is  200  length 
units,  and  the  hollow  cone  antenna  pattern  is  used.  The  signature 
sampling  increment  along  the  trajectory  is  3  length  units,  and  each 
sub-area  of  the  sea  surface  has  a  width  of  6  length  units.  The 
number  of  sub-areas  is  21  by  13  with  21  sub-areas  located  transverse 
to  the  trajectory  direction.  With  pattern  cone  half-angle  of  60°, 
the  three  different  signatures  in  Figures  66  through  68  correspond 
to  three  different  antenna  approach  angles.  The  signatures  in  Fig¬ 
ures  68  through  71  represent  the  signatures  for  pattern  cone  half 
angle  of  30°.  The  starting  location  fSRF)  of  the  radar  from  the 
center  of  the  ship  has  to  vary  drastically  from  one  antenna  approach 


ang Ip  to  another.  This  is  caused  hv  the  substantial  size  difference 
between  thp  top  and  bottom  cylinders.  For  pxamplp,  when  thp  antpnna 
has  a  broadside  incidpnce  ISAH=0°1,  the  ship-scattpred  and  the  ship- 
sea  double-bounce  signatures  are  very  close  to  each  other  in  space; 
however,  when  the  antenna  approaches  head-on  (SAH=90°1,  the  above 
two  signatures  are  very  far  apart.  Thus,  in  order  to  show  both 
signatures  on  a  single  run,  both  the  antenna  starting  location  and 
the  distance  it  travels  need  to  vary  depending  on  the  antenna  ap¬ 
proach  angle.  It  can  be  seen  that  the  ship  scattered  signatures 
are  easily  detectable,  independent  of  the  antenna  approach  angle 
or  the  pattern  cone  half  angle.  The  ship-sea  double-bounce  signa¬ 
ture,  however,  is  only  detectable  when  the  antenna  approaches  from 
broadside.  This  characteristic  so  far  has  not  varied  for  any  of 
the  models  that  have  been  used  in  this  report.  It  should  be  pointed 
out  that  the  scattered  signatures  due  to  the  ship  alone  is  believed 
to  be  a  very  crude  approximation  to  the  real  scattered  signature. 

This  is  based  on  the  fact  that  there  are  many  substructures  on  the 
ship  that  can  scatter  energy.  Figure  7?  demonstrates  that,  to  im¬ 
prove  the  simplified  U.  S.  S.  Missouri  Model,  several  objects  such 
as  blocks,  cylinders  and  plates  should  be  added  onto  the  decks  of 
the  ship.  The  GT(1  technique  can  still  he  effectively  applied  to 
find  the  scattering  from  such  a  complicated  ship  structure  at.  the 
cost  of  some  extra  complexity  in  the  computer  program. 

In  actuality,  the  ship  contribution  to  the  signature  from  either 
the  more  complex  model,  or  the  real  ship,  will  be  considerably  dif¬ 
ferent  in  duration,  magnitude,  and  structure  from  those  shown  In 
the  nrecedinq  fiqures.  Our  modelling  of  the  hull  and  sea  surface 
is  a  reasonably  good  approximation  to  the  real  situation,  and  it 
may  be  expected  that  the  clutter  and  double  bounce  signatures  will 
he  quite  similar  to  those  obtained  here.  It  is  most  unlikely  how¬ 
ever  that  a  real  ship  would  exhibit  the  long  expanse  of  unbroken 
deck  and  cabin  wall  needed  to  form  the  continuous  corner  reflector 
which  is  the  major  contributor  to  a  ship  signature.  Instead,  the 
corner  would  be  broken  up  by  obstructions  and  augmented  by  reflec¬ 
tions  from  other  structures,  producing  an  irregular  signature  of 
somewhat  longer  duration  (especially  for  a  bow  or  stern  approach). 

It  is  probable,  however,  that  the  real  ship  signature  would  not 
ever  be  significantly  larger  in  magnitude  than  the  largest  signa¬ 
tures  shown  here,  because  single,  double  or  triple  bounce  reflected 
fields  are  all  essentially  image  fields.  In  any  event,  no  attempt 
should  be  made  to  evaluate  system  performance  on  the  basis  of  the 
signatures  in  this  report,  because  of  the  misleading  simplicity 
of  the  ship-alone  component. 
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r  i quro  5.1.  Tot.il  scattered  signature  from  the  sea  and  a  ship  with 
conical  hull.  Hollow  cone  antenna  pattern  is  used.  For  the 
geometrical  configuration  of  Figure  45,  SAH-n  .  Antenna 
pattern  cone  half  angle  ATA~RP  .  Ship's  center  is 
located  at  RF^lOO  length  units  with  RF  being  the 
distance  traveled  hv  the  antenna  from  its  starting 
location.  Ship's  dimensions  are  A1-5,  A?-*, 

R?=-vv  7fY-'1d ,  TA,  "TY'^'n  length  units. 
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Figure  S4 .  Total  scattered  siqnature  from  the  sea  and  a  ship  with 
conical  hull.  Hollow  cone  antenna  pattern  is  used.  SAH«4F°, 
ATA=60°.  Other  dimensions  are  the  same  as  those 
shown  in  Figure  F3. 
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Figure  55.  Total  scattered  signature  from  the  sea  and  a  ship 
with  conical  hull.  Hollow  cone  antenna  pattern  is  used. 
SAH=90°,  ATA=60°.  Other  dimensions  are  the  same  as 
those  shown  in  Figure  53. 


Figure  56.  Total  scattered  signature  for  the  same  geometrical 
configuration  as  that  for  Figure  53.  Here  the  antenna  pattern 
cone  half  angle  ATA»30°.  Antenna  approach  angle  SAH»0  . 
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Figure  57.  Total  scattered  signature  from  the  sea  and  a  ship 
with  conical  hull.  Hollow  cone  antenna  pattern  is  used. 
SAH=45  ,  ATA=30°.  Other  dimensions  are  the  same  as 
those  shown  in  Figure  53. 
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Fiqure  F8.  Total  scattered  siqnature  from  the  sea  and  a  ship 
with  conical  hull.  Hollow  cone  antenna  pattern  is  used. 
SAH=«0  ,  ATAa30°.  Other  dimensions  are  the  same  as 
those  shown  in  Fiqure  5.1. 
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Figure  59.  Total  scattered  signature  from  the  sea  and  a  ship 
with  conical  hull.  Hollow  cone  antenna  pattern  is  used. 
SAH=0°,  ATA=60°,  polarization  is  horizontal. 

Other  dimensions  are  the  same  as  those 
shown  in  Figure  53. 
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Fiqure  FO.  Total  scattered  siqnature  from  the  sea  and  a  ship 
with  conical  hull.  Hollow  cone  antenna  pattern  is  used. 
SAH=4F°,  ATA=AO  ,  polarization  is  horizontal. 

Other  dimensions  are  the  same  as  those 
shown  in  Fiqure  S3. 
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f iqure  ai  .  Total  scattered  siqnature  from  the  sea  and  a  ship  with 
conical  hull.  Hollow  cone  antenna  pattern  is  used. 

SAH=°rr,  ATA=fiO  ,  polarization  is  horizontal. 

Other  dimensions  are  the  same  as  those 
shown  in  Figure  FT. 
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re  62.  Total  scattered  signature  from  the  sea  and  a  ship  with 
conical  hull.  Hollow  cone  antenna  pattern  is  used. 

SAH=0°,  ATA=30°,  polarization  is  horizontal. 

Other  dimensions  are  the  same  as  those 
shown  in  Figure  53. 
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Figure  63.  Total  scattered  signature  from  the  sea  and  a  ship  with 
conical  hull.  Hollow  cone  antenna  pattern  is  used. 

SAH=45  ,  ATA=30  ,  polarization  is  horizontal. 

Other  dimensions  are  the  same  as  those 
shown  in  Figure  63. 
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Figure  66.  Ship-sea  scattered  signature  from  the  simplified 
U.S.S.  Missouri  model.  Pattern  cone  half  angle  ATA=60  . 
Antenna  approach  angle  SAH=0°.  RF  is  the  distance  be¬ 
tween  the  antenna  and  its  starting  position. 

Vertical  scale  represents  the  received 
voltage  in  decibles.  Ship’s  center  is 
located  at  RF=450  length  units  and 
ship's  dimensions  are  shown  in 
Figure  65. 
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Fiqure  fi7.  Ship-sea  scattered  signature  from  the  simplified 
U.S.S.  Missouri  model.  SAH=4*°,  ATA=fiO  .  RF  is  the 
distance  between  the  antenna  and  its  startino  position  and 
ship's  center  is  located  at  RF=^4n  length  units.  Ship's 
dimensions  are  shown  in  Figure  A*’. 
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Figure  69.  Ship  and  sea  scattered  signature  from  the  simplified 
U.S.S.  Missouri  model.  Pattern  cone  half  angle  ATA=30  , 

antenna  approach  angle  SAH=0  .  RF  is  the  distance  .1 

between  antenna  and  its  starting  position.  Ship's 

center  is  located  at  RF=775  length  units.  Ship's  ■ 

dimensions  are  shown  in  Figure  65.  i 
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Mqure  7o.  Ship-sea  scattered  siqnature  from  the  simplified 
U.S.S.  Missouri  model.  SAH*AS,  ATA=S0  .  RF  is  the 
distance  between  the  antenna  and  its  startinn  position, 
and  ship's  center  is  located  at  RF=i?nn  length  units. 
Ship's  dimensions  are  shown  in  Figure  FF. 
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Figure  71.  Ship-sea  scattered  signature  from  the  simplified 
U.S.  Missouri  model.  SA H»Qf>  ,  ATA*3f>  .  RF  is  the 
distance  between  the  antenna  and  its  starting  position, 
and  ship's  center  is  located  at  RF»1?00  length  units. 
Ship's  dimensions  are  shown  in  Figure  65. 


CHAPTER  VII 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  objective  of  this  research  has  been  to  develop  a  simple 
model  for  the  interaction  hetween  a  radar  system  and  a  ship  at  sea 
which  is  adaptable  to  short  ranqe  encounters,  and  which  can  be  easily 
elaborated  to  handle  fairly  complex  ship  structures.  In  the  most 
qeneral  form  considered  in  this  report,  the  superstructure  of  the 
ship  is  modelled  by  an  elliptic  cylinder,  and  the  hull  is  modelled 
by  a  portion  of  an  elliptic  cone.  The  sea  is  modelled  by  an  array 
of  rectangular  patches,  each  characterized  by  a  polarization  depen¬ 
dent  bistatic  radar  cross-section.  The  radar  system  response  is 
represented  by  an  encounter  signature,  i.e.,  by  the  detector  voltage 
as  a  function  of  the  radar  position.  The  total  signature  is  then 
found  from  the  (phasorl  sum  of  the  contributions  from  each  of  the 
dominant  scattering  mechanisms.  These  are  the  clutter  (back-scat¬ 
tered)  signal  from  the  sea  alone;  the  reflected  and  diffracted  fields 
from  the  ship  structure  alone;  and  the  double-bounce  contribution 
produced  by  the  fields  scattered  by  the  ocean  surface  onto  the  ship 
structure  and  then  back  to  the  radar  system  (or  vice  versa). 

The  fields  scattered  by  the  various  component  parts  of  the 
ship  are  calculated  by  the  Geometrical  Theory  of  Diffraction  (GTD), 
combined  with  a  very  fast  numerical  search  technique  which  was  de¬ 
veloped  to  track  the  reflection  and  diffraction  points  on  the  ship’s 
structure. 

From  the  perspective  afforded  by  the  detailed  study  of  the 
computational  procedures  described  in  the  preceding  chapters,  the 
advantages  of  this  particular  ship-sea  model  should  be  readily  ap¬ 
parent.  In  the  first  place,  since  the  GTD  assigns  a  definite  ray 
path  to  each  contribution  to  the  total  signature,  it  is  relatively 
simple  to  adapt  the  model  to  a  variety  of  features  found  in  oper¬ 
ating  radar  systems,  e.g.,  ranqe  gates,  doppler  filters,  pulse-to- 
pulse  frequency  aqilit.v,  etc.  Secondly,  because  the  total  scattered 
field  from  each  patch  is  radiated  from  the  center  of  the  patch, 
the  location  of  which  is  known,  only  one  point  f that  on  the  ship's 
surface)  must  be  located  to  define  the  ray  path.  Since  a  large 
number  of  paths  must  be  computed  (one  for  each  patch  at  each  position 
of  the  radar),  this  feature  of  the  model  greatly  reduces  the  overall 
computation  time  compared  to  any  model  in  which  two  unknown  reflec¬ 
tion  points  per  path  must  be  tracked.  Thirdly,  since  the  actual 
bistatic  scattering  cross-section  of  the  ocean  is  used  to  determine 
the  strength  of  the  field  scattered  by  each  patch,  both  the  clutter 
signature  and  the  double  bounce  signature  should  have  the  same  stat¬ 
istical  properties  (average  power,  correlation  function,  spectrum. 


lOfi 


first  probability  distribution)  as  do  the  actual  signatures.  Thus, 
one  may  use  the  model  signatures  to  evaluate,  for  example,  proposed 
antenna  patterns,  detection  characteristics,  transmitted  signal 
designs,  polarization  diversity  schemes,  etc.  For  example,  even 
though  the  random  clutter  signal  may  occasionally  attain  large  values, 
i.e.,  in  relation  to  the  ship  return,  it  could  be  distinguished 
from  the  ship  return  by  the  fact  that  it  can  sustain  large  values 
only  for  very  short  durations. 

Although  only  relatively  few  signatures  were  actually  calcu¬ 
lated  in  preparing  this  report,  one  can  already  observe  a  variety 
of  signature  characteristics  which  may  be  of  value  in  future  system 
design.  For  example,  it  is  clear  that  the  "double-bounce"  contri¬ 
bution  may  be  significantly  larger  than  the  clutter,  particularly 
at  broadside  and  bow  or  stern  aspects,  and  that  the  early  position 
of  its  maximum  response  is  determined  by  the  waterline  and  hull 
rather  than  the  superstructure  details.  On  the  other  hand,  (for 
pencil  beam  antennas)  there  are  aspects  near  45°  for  which  neither 
the  double-bounce  nor  the  ship  signature  is  detectable  in  the  clut¬ 
ter.  When  a  hollow  cone  antenna  pattern  was  used,  the  presence 
of  the  ship  was  always  detectable  even  though  the  clutter  level 
was  higher.  However,  in  this  case  one  could  observe  either  a  single 
peak  or  a  double-peak  in  the  overall  signature  depending  again  on 
the  aspect  angle  of  the  trajectory.  As  mentioned  on  nage  85,  how¬ 
ever,  one  should  not  use  the  signatures  of  the  type  developed  in  this 
report  to  evaluate  system  performance,  because  of  the  extreme  sim¬ 
plicity  of  modelling  the  superstructure. 

On  the  basis  of  these  brief  observations,  and  the  other 
results  and  the  signatures  of  Chapter  VII,  it  may  be  concluded  that 
the  simple  ship-sea  scattering  model  developed  here  offers  a  con¬ 
venient  method  for  estimating  the  performance  of  a  radar  system 
traversing  a  ship  at  sea.  It  should  be  useful  to  the  system  de¬ 
signer  in  evaluating  the  effects  of  such  system  parameters  as  an¬ 
tenna  pattern,  polarization,  radar  trajectory,  range  gate,  detection 
criteria,  etc.,  when  the  modelling  of  the  superstructure  has  been 
augmented  by  appropriate  forms  (plates,  cones,  cylinders,  blocks, 
etc.)  all  of  which  can  be  handled  by  the  technique  of  the  GTD  im¬ 
plemented  here.  In  addition,  the  sea  surface  model  can  be  elab¬ 
orated  to  include  non-random  components  (e.g.,  the  bow-wave  struc¬ 
ture  always  associated  with  ships  in  motion,  which  may  be  the  domi¬ 
nant  surface  disturbance  in  smooth  water).  Alternatively,  one  could 
use  measured  signatures  from  an  isolated  model  of  the  ship,  in  con¬ 
junction  with  computed  signatures  for  the  clutter  and  double  bounce 
contributions  to  obtain  an  overall  signature. 
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appendix  a 

PENCIL  BEAM  RADIATION  FROM  A  RECTANGULAR  APERTURE 


An  antenna  with  the  rectanqular  aperture  and  coordinate  sys¬ 
tem  show^in  Fiqure  emits  a  radiated  far  2one  electric  field 
qiven  hy 


L/P  W /P 

/  J  [sine  J  fx’.z'l-cosecos*  J  f x '  z ' >1 
-L IP  -W/2  z  x  *  J 


e jk  f  x '  $  i  n0cos*+z '  cose )  dx<  ^  , 
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UP  W/2 

f  J 

-L/2  -W/2 


j  (x.<zi^jk(x,sinecos*+2'cose) 
x  *  dx ' dz ' 


(60) 

when  the  origin  of  the  coordinates  is  taken  at  the  center  of  the 
aperture,  and  the  spherical  coordinate  system  for  the  radiated  field 
is  as  shown. 


z 


Figure  73.  Source  coordinate  system. 
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Since  the  assumed  antenna  current  distribution  is  Gaussian, 
the  z-directed  aperture  current  is  given  by 


.Hx',*') 
Jx  =0 


V 


-n(x'/aV 


■n(z*/aV 


(61) 


where  a  is  a  constant  parameter  that  determines  the  beamwidth  and 
has  units  of  length.  J  is  a  reference  surface  current  density  in 
units  of  amperes  per  meter.  Since  the  Gaussian  source  function  is 
diminishingly  small  outside  the  limits  of  integration,  the  inte¬ 
grals  can  be  approximated  by  letting  the  limits  go  to  infinity, 

with  the  result  that  „  .?  .  ,? 

(ka  cos 0)  (ka  sinflcosj)! 

=  J  a2  e-Jkr  sine  <  * 


"0  4  nr 


•  j  J  sin6 

4nr  0 


0+sin^ecos?4i) 


(62) 


=  0. 


Similarly,  for  an  x-directed  Gaussian  aperture  current,  the  elec 
trie  fields  are  given  by 


zp H  j_a?  e"jkr 
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cosftcos4>  e 
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(cos2e+sin2ecos2<i») 
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(64) 


Since  the  antenna  generates  a  narrow  broadside  (ft=Q0°,  <b=°f)°l  beam, 
and  cosfi  cos<h  is  vanishingly  small  near  the  broadside  direction, 
is  negligible  for  the  x  polarized  source  distribution. 

To  conclude,  the  fields  that  are  incident  upon  the  sea  surface 
and  the  ship  are  given  by  Equations  (69)  and  (64)  for  a  vertically 
(J  )  and  horizontally  (J  )  polarized  antenna,  respectively.  In 
the  report,  all  the  above  calculations  have  assumed  0  a  =1  ampere¬ 
meter. 
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APPFNWX  B 


BISTATIC  SCATTERING  MATRIX  ELEMENTS 
(slightly  rough  surface) 


The  histatic  scattering  matrix  elements  for  a  slightly  rough 
surface  with  homogeneous  dielectric  material  are  given  by1 
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In  the  above  expressions,  e  ,  the  relative  permittivity,  may  be 
either  real  or  complex  (indicating  lossless  or  lossy  material), 
and  may  be  either  greater,  equal  to,  or  less  than  unity  in  magni¬ 
tude.  The  angles  o  0  and  <J>  are  defined  in  Figure  8.  The  first 
index  n  of  onr  is  the  scattered  polarization  state  and  the  second 
index  c  is  tne  Incident  polarization  state.  For  a  perfect  conduct¬ 
ing  surface 


110 


APPENDIX  C 

BISTATIC  SPECULAR  REFLECTION  COEFFICIENTS 
(very  rough  surface) 


The  specular  reflection-coefficients  for  a  surface  with  larqe 
scale  roughness  are  given  by71 
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The  quantities  R„(v)  and  RA(v)  are  the  Fresnel  reflection  coeffl 
dents  defined  as 
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cosv 


(67b) 


with  A.,  A  and  <|>  defined  in  Fiqure  A.  Note  that  special 
care  i£  required  in  evaluating  these  coefficients  in  the  back-scat- 
terinq  direction.  For  this  reason,  a  somewhat  different  procedure 
for  evaluating  these  coefficients  are  used  in  the  computer  programs 
in  this  report. 
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APPENDIX  D 

COMPARISON  OF  THE  MEASURED  AND  CALCULATED  FIELDS  SCATTERED 

BY  A  CIRCULAR  CONE 


In  order  to  validate  the  calculation  for  the  reflected  field 
from  an  elliptic  cone,  the  most  appropriate  verification  Is  to  com¬ 
pare  the  results  with  experimental  measurements,  since  neither  cal¬ 
culated  nor  measured  results  are  available  In  the  standard  literature. 
However,  because  the  elliptic  cone  is  very  difficult  to  fabricate, 
the  experimental  data  reported  here,  and  consequently  the  theoreti¬ 
cal  solution  derived  In  this  appendix.  Is  restricted  to  the  case 
of  the  circular  cone.  The  experimental  results  were  obtained  bv 
Dr.  Fdward  Pelton  In  the  ESL  (ElectroSclence  Laboratory!  microwave 
anechoic  chamher  usinq  S-hand  setups.  The  field  scattered  by  the 
circular  cone  Is  received  by  a  rectanqular  horn  antenna  and  detected 
by  a  Scientific  Atlanta  receiver. 

The  experimental  conf iquration  of  the  circular  cone  is  illus¬ 
trated  in  Figure  74.  a  half-wave  dipole  which  is  located  in  the 
near  zone  of  the  cone  illuminates  the  cone  with  either  vertical 
or  horizontal  polarization,  and  the  near-zone  field  is  measured 
azimuthally  around  the  cone.  Since  the  measured  result  Includes 
not  only  the  reflected  field  but  also  the  incident  field  plus  the 
surface  diffracted  creeping  wave.  It  is  essential  to  calculate  all 
these  fields  in  order  to  have  an  adequate  comparison.  The  reflected 
field  in  the  deep  lit  region  (away  frun  the  shadow  boundary)  is 
calculated  by  the  same  Geometrical  Optics  technique  that  Is  employed 
for  calculating  the  field  reflected  from  the  ship's  hull.  The  re¬ 
flected  field  near  the  shadow  boundary  (transition  region)  and  the 
surface  diffracted  creeping  wave  are  both  calculated  by  the  uniform 
GTD  method  using  the  formulation  of  Pathalc' . 

Refore  calculating  the  creeping  wave  field  it  is  essential 
to  accurately  locate  the  geodesic  ray  path  for  each  creeping  wave 
on  the  conical  surface,  the  geodesic  curve  on  a  surface  is  the 
shortest  path  between  two  points.  For  given  locations  of  the  source 
'x  y  z  )  and  observation  point  fx  ,y  ,z  ),  the  qeodesic  ray  path 
onsa  circular  cone  can  he  found  as  follows: 

Let  the  incident  and  the  outgoing  vector  tanqent  points  be  denoted 
by  01  and  0,,  respectively,  as  illustrated  in  Figure  7R,  with 


Q,  »  P| cos  V| x  +  Pj  sin  v^y  +  p]CotPo7 
*  P?cos  v?x  ♦  p2sin  v?y  +  p2cotflQz 


PERSPECTIVE  VIEW  : 
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Figure  74.  Geometry  and  scattering  mechanisms 
for  the  circular  cone. 


PERSPECTIVE  VIEW  TOP  VIEW 


Figure  75.  Diffraction  mechanism  for  creeping  wave. 


where  p. ,  p?,  v.,  v?  are  the  unknowns  to  be  solved  for,  with  p, 
heing  horizontal  ranil  of  the  cone  and  v^  being  angular  positions 
measured  from  the  x-axis;  0Q  is  the  conenalf  angle. 

The  incident  vector  f  and  the  outgoing  vector  R  can  be  described 
by 

^  A  A  ^ 

I  *  (p|Cos  v1-xslx  +  fp|Sin  V|-.vsly  +  (pj  cotflQ-zslz 

v  A  A  A 

R  ■  fx0-p?cos  v?)x  ♦  (,y0-PjSln  v?)y  +  (zo-P,>cot0o)z  .  (70) 

The  normal  to  the  surface  Is  given  by  (see  Equation  (46)) 

A  A  A 

«  cos  v^  coteQx  +  sin  v^  cot  BQy  -  z.  (71) 

Since  the  ray  path  has  to  be  everywhere  perpendicular  to  the  surface 
normal,  at  the  Incident  point  one  must  satisfy 

I  •  Nj  *  0 

which  results  in  the  equation 


At  the  outgoing  tangent  point  one  should  have 
R  •  N?  »  0 

which  yields  the  equation 

zQ  -  cot  e0(xQ  cos  V/j  +  ys  sin  v?)  »  0.  (73) 

Equations  (72)  and  (73)  allow  one  to  solve  for  the  angular  positions 
v,  and  v?  independent  of  the  distances  p,  and  p„,  which  can  be  de¬ 
terminedly  the  following  method.  1  £ 

If  a  circular  cone  Is  cut  open  along  a  generating  line  and 
is  un-rolled  onto  a  flat  surface,  one  obtains  the  so-called  developed 
surface,  which  is  illustrated  in  Figure  76.  The  distances  from 
the  tip  of  the  cone  to  the  source  and  observation  points  are  given 
by 


Figure  76.  Developed  cone  surface. 


and 

f~2  ?  V 

h=J*  o  +  *o+2o’ 

respectively,  with  unit  direction  vectors  given  by 

*1  ■  <v  ♦  v  *  V>/li 

A  A  A  A 

*2  *  (V  +  v  +  zoz)/l2 

The  generating  lines  for  the  cone  at  the  two  tangent  points  are 
described  by  the  following  unit  vectors 

-A  A  /A  /A 

r^  =  cos  Vj  sineQx  +  sin  v^  s1n60y  +  cos8Qz 

/A  A  A  A 

r?  »  cos  v?  sin60x  +  sin  v^  s1neQy  +  cos0Qz. 

The  angles  8^,  as  shown  in  Figure  76,  are  given  by 
61  =  (v?-Vj)sinfto 

A  /A 

B-j  =  arccosf  «r|l 

A  A 

R?  =  arccosi£,?*r?l 

Thus,  the  ray-path  distance  between  the  source  and  observation  points 
is 


s  =  +  4 ' 2%lh.  cos{81+B2+ft1)]55,  f 74) 

and  the  angles  a,  and  c u,  as  illustrated  in  the  same  figure,  are 
given  by 


aj  a  arccos[(s^+s2-«|)/2a1s] 
c?  a  *“«i  -(e1+e2+e1) 

Finally,  pj  and  p^  can  be  obtained  from 


sine 


*  sTnTo[+8j') 


sin  0. 


(75) 


P? 


*?  s1n<^> 
sinf  0,+Bpl 
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After  obtaining  the  locations  of  the  two  tanqent  points  (P.,v,) 
anil  (p->*  v.>) ,  one  is  a  We  to  solve  for  the  constants  c.  and1 c-1 in 
the  geodesic  equation  ,  1 

Pi  -  c1sec[(vi-c?)sin<>0]  .  (76) 


The  surface  diffracted  field  of  the  creeping  wave  (see  Figure 
y7)  has  been  derived  by  Pathak  and  is  given  by''' 


C  JkS 


s(oc+sl 


(77) 


'  '  K  A  A  AAA 

where  and  b.,  are  the  binormal  unit  vectors  (b.-t.xn. ,  b?«t?xn,,l , 
and  F  *0^1  isthe  arbitrary  polarized  incident  field  at  point  ; 

Jdn^/dr^  ’  is  the  rav  spreading  factor; 

0  is  the  radius  of  the  spherical  wave  front  at  Q.,  and  is  qiven 
r  by  oc*s’*t  with  t  being  the  geodesic  path  lenoth. 


A  A 

b,  n, 


Figure  77.  Configuration  of  surface  diffracted 
creeping  wave. 


The  diffraction  coefficient  D  and  D.  for  the  soft  and  hard  boundary 
conditions,  respectively,  aresqiven  By 

f  i  /a  > 

Ds  =  e'jkt [l-F(xs1]  +  p$^S1  '  f7R) 

h  h 


where 


m(0. )  »  kYQ<) 


»,(«,)  '  -J—1 - ?- 

y  C|COS0Q  cos  w. 


where  wi  =  ^  -  0^  -  8^ 


is  the  equation  for  the  radius  of  curvature  of  the  geodesic  ray 
path  and  is  obtained  in  a  fashion  similar  to  Equation  (53).  The 
constant  c.  Is  the  same  constant  as  that  shown  in  Equation  (76). 
The  distance  parameter  ($s)  from  Qj  to  Q2  is  defined  by 


c  ??  m(t') 

£*•  n  ojT) 


dt  *  . 


For  a  circular  cone,  E  is 


Es  =  ^  kp1  sinQ^  |  Vp-v^  |cos9/,^Pc 


where 


nl  =  arcsin  |s-"sTn60  sin|(vv1)s1ne0||. 

The  function  F(xs)  Is  formulated  In  terms  of  the  Fresnel  integral, 
F(xs)  >  2jJ?  e'jx  J  e'jT'  dx  (80) 


where 


s  _  ks's 


s'+s  ?m(Q1)m(Q2) 


The  function  P(£ )  is  the  Pekeris*  caret  function  and  is  identical 
s 

^  24  23 

to  the  definitions  of  Logan  and  Pathak  . 

The  qeometrical  optics  field  will  not  blend  accurately  with 
the  field  in  the  shadow  region  unless  the  geometrical  reflected 
field  is  modified  in  the  transition  region  according  to  (see  Fig¬ 
ure  781 

rr(pl  *  1  -fe—— —  e"jks  (81) 


(p'+sU^+s) 


.  1  1  ,  1 
where  =_  +_ 


with  pj  being  the  principal  radii  of  curvature  of  the  reflected 

wavefront  at  the  reflection  point  Q„  for  a  plane  wave  incidence, 
and  s'  being  the  distance  between  tne  source  point  and  QR 


Figure  78.  Configuration  of  curved  surface  reflection. 


The  dyadic  reflection  coefficient  is  given  by 


R  *  R  eAex+  Rh  ej 


1 


where 


-4  e 


RS"-'- 
h 


•JCf.1)3/!?/  e~jTr/4 

\?-fn  K 


1-F(x*')+  P$f^) 


) 


with  xl  =  2  cosV 


ll  =  -?m  h1/3  cose1 


P  tfv  F  „ .  r 

^  VVcosV  1 

nr—  J  •  h  1 

J 

Here  p  (Q„)  is  the  radius  of  curvature  of  the  surface  at  Q„  in  the 
plane  of  incidence;  whereas,  p.  (QR)  is  the  radius  of  the  surface 
in  the  orthogonal  plane. 


The  reflected  field  from  the  circular  cone  is  calculated  by 
Equation  (81)  if  the  field  point  is  close  to  the  shadow  boundary 
or  by  Equation  (1)  if  the  field  point  is  in  the  deep  lit  region. 

The  total  field  is  then  composed  of  the  incident  field  and  the  re¬ 
flected  field  plus  the  surface  diffracted  creeping  wave  calculated 
by  Equation  (77).  The  comparison  between  this  calculated  total  field 
and  the  measured  result  from  the  geometry  given  in  Figure  74  is 
shown  in  Figure  79.  It  can  be  seen  that  excellent  agreement  is 
obtained  between  the  GTD  analysis  and  the  measurement.  Consequently, 
one  can  conclude  that  the  reflected  field  from  the  elliptic-cone 
shaped  ship  hull  has  been  correctly  calculated. 
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FIELD  (  OB)  ,  FIELD  (db  ) 


♦  -ANGLE  tOEGREES) 


RF  .  32.32X 
ZF  •  I8.I5X 

LENGTH  OF  CONE  &  =  14.56  X 


Figure  79.  Field  scattered  from  the  circular  cone.  (Solid 
curve  is  the  calculation,  dashed  curve 
is  the  measurement. 
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